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Abstract 

We report on work to increase the number of well-measured Type la supernovae (SNe la) at high 
redshifts. Light curves, including high signal-to-noise HST data, and spectra of six SNe la that were 
discovered during 2001 are presented. Additionally, for the two SNe with z > 1, we present ground- 
based J-band photometry from Gemini and the VLT. These are among the most distant SNe la for 
which ground based near-IR observations have been obtained. We add these six SNe la together 
with other data sets that have recently become available in the literature to the Union compilation 
(Kowalski et al. 2008). We have made a number of refinements to the Union analysis chain, the 
most important ones being the refitting of all light curves with the SALT2 fitter and an improved 
handling of systematic errors. We call this new compilation, consisting of 557 supernovae, the Union2 
compilation. The fiat concordance ACDM model remains an excellent fit to the Union2 data with the 
best fit constant equation of state parameter w = — 0.997;l^Q 054(stat)^g Qg2(stat -|- sys together) for a 
flat universe, or u; = — 1.035lQ;Q5g(stat)^Q Qg7(stat -I- sys together) with curvature. We also present 
improved constraints on w{z). While no significant change in w with redshift is detected, there is still 
considerable room for evolution in w. The strength of the constraints depend strongly on redshift. In 
particular, at z > 1, the existence and nature of dark energy are only weakly constrained by the data. 
Subject headings: Supernovae: general — cosmology: observations — cosmological parameters 
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1. INTRODUCTION 

Type la supernovae (SNe la) are an excellent tool for 
probing the expansion history of the Universe. About a 
decade ago, combined observations of nearby and distant 
SNe la led to the discovery of the accelerating universe 
(Perlmutter et al. 1998; Garnavich et al. 1998; Schmidt 
et al. 1998; Riess et al. 1998; Perlmutter et al. 1999). 

Following these pioneering efforts, the combined work 
of several different teams during the past decade has pro- 
vided an impressive increase in both the total number 
of SNe la and the quality of the individual measure- 
ments. At the high redshift end (z > 1), the Hubble 
Space Telescope (HST) has played a key role. It has suc- 
cessfully been used for high-precision optical and infrared 
follow-up of SNe discovered from the ground (Knop et al. 
2003; Tonry et al. 2003; Barris et al. 2004; Nobili et al. 

2009) , and, by using the Advanced Camera for Surveys 
(ACS), to carry out both search and follow-up from space 
(Riess et al. 2004, 2007; Kuznetsova et al. 2008; Daw- 
son et al. 2009). At the same time, several large-scale 
ground-based projects have been populating the Hubble 
Diagram at lower redshifts. The Katzman Automatic 
Imaging Telescope (Filippenko et al. 2001), the Nearby 
Supernova Factory (Copin et al. 2006), the Center for As- 
trophysics SN group (Hicken et al. 2009a), the Carnegie 
Supernova Project (Hamuy et al. 2006; Folatelli et al. 

2010) , and the Palomar Transient Factory (Law et al. 
2009) are conducting searches and/or follow-up for SNe 
at low redshifts (z < 0.1). The SN Legacy Survey (SNLS) 
(Astier et al. 2006) and ESSENCE (Miknaitis et al. 2007; 
Wood-Vasey et al. 2007) are building SN samples over the 
redshift interval 0.3 < 2; < 1.0, and the SDSS SN Survey 
(Holtzman et al. 2008; Kessler et al. 2009) is building 
a SN sample over the redshift interval 0.1 < z < 0.3, 
a redshift interval that has been relatively neglected in 
the past. These projects have discovered ~ 700 well- 
measured SNe. The number of well-measured SNe be- 
yond z ~ 1 is approximately 20 and is comparatively 
small. 

Kowalski et al. (2008) (hereafter K08) provided a 
framework to analyze these and future datasets in a 
homogeneous manner and created a compilation, called 
the "Union" SNe la compilation, of what was then the 
world's SN data sets. Recently, Hicken et al. (2009b) 
(hereafter H09) added a significant number of nearby 
SNe to a subset of the "Union" set to create a new 
compilation, and similarly the SDSS SN survey (Kessler 
et al. 2009) (hereafter KS09) carried out an analysis of a 
compilation including their large intermediate- 2: data set 
(Holtzman et al. 2008). When combined with baryon 
acoustic oscillations (Eisenstcin et al. 2005), the H09 
compilation leads to an estimate of the equation of state 
parameter that is consistent with a cosmological constant 
while KS09 get significantly difi'erent results depending 
on which light curve fitter they use. 

An important role for SNe la beyond 2; ~ 1, in addition 
to constraining the time evolution of w, is their power to 
constrain astrophysical effects that would systematically 
bias cosmological fits. Most evolutionary effects are ex- 
pected to monotonically change with redshift and are not 
expected to mimic dark energy over the entire redshift 
interval over which SNe la can be observed. Evolution- 
ary effects might also have additional detectable conse- 



quences, such as a shift in the average color of SNe la 
or a change in the intrinsic dispersion about the best fit 
cosmology. 

Interestingly, the most distant SNe la in the Union 
compilation (defined here as SNe la with z>l.l) are al- 
most all redder than the average color of SNe la over the 
redshift interval 0.3 < z < 1.1. The result is unexpected 
as bluer SNe la at lower redshifts are also brighter (Tripp 
1998; Guy et al. 2005) and should therefore be easier to 
detect at higher redshifts. Possible explanations for the 
redder than average colors of very distant SNe la range 
from the technical, such as an incomplete understanding 
of the calibration of the instruments used for obtaining 
the high redshift data, to the more astrophysically inter- 
esting, such as a real lack of bluer SNe at high redshifts. 

The underlying assumption in using SNe la in cos- 
mology is that the luminosity of both near and distant 
events can be standardized with the same luminosity 
versus color and luminosity versus light curve shape re- 
lationships. While drifts in SN la populations are ex- 
pected from a combination of the preferential discovery 
of brighter SNe la and changes in the mix of galaxy types 
with redshift (Howell et al. 2007) — eff'ects that will aff'ect 
different surveys by differing amounts a lack of evo- 
lution in these relationships with redshift has not been 
convincingly demonstrated given the precision of current 
data sets. This assumption needs to be continuously ex- 
amined as larger and more precise SN la data sets become 
available. 

In this paper, we report on work to increase the num- 
ber of well-measured distant SNe la by presenting SNe la 
that were discovered in ground based searches during 
2001 and then followed with WFPC2 on HST. Two of 
the new SNe la are at z ~ 1.1 and have high-quality 
ground-based infrared observations that were obtained 
with ISAAC on the VLT and NIRI on Gemini. This 
paper is the first paper in a series of papers that will 
provide a comparable sample of 2; > 1 SNe la to the SNe 
now available in the literature. The SNe la in this se- 
ries of papers were discovered in 2001 (this paper), 2002 
(Suzuki et al. in preparation) and from 2005 to 2006 
during the Supernova Cosmology Project (SCP) cluster 
survey (Dawson et al. 2009). 

The paper is organized as follows. In Section 2, we 
describe the SN search and the spectroscopic confirma- 
tion, while Sections 3, 4 and 5 contain a description of 
the follow-up imaging and the SN photometry. The light 
curve fitting is described in Section 6. In Section 7 we 
update the K08 analysis both by adding new data and 
by improving the analysis chain. The paper ends with a 
discussion and a summary. 

2. SEARCH, DISCOVERY AND SPECTROSCOPIC 
CONFIRMATION 

The SNe were discovered during two separate high- 
redshift SN search campaigns that were conducted dur- 
ing the Northern Spring of 2001. The first campaign 
(hereafter Spring 2001) consisted of searches with the 
CFH12k (Cuillandre et al. 2000) camera on the 3.4 m 
Canada-Francc-Hawaii Telescope (CFHT) and the MO- 
SAIC II (MuUer et al. 1998) camera on the 4.0 m Cerro 
Tololo Inter-American Observatory (CTIO) Blanco tele- 
scope. The second campaign (hereafter Subaru 2001) 
was done with SuprimeCam (Miyazaki et al. 2002) on 
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the 8.2 m Subaru telescope. AU searches were "classical" 
searches (Perlniutter et al. 1995, 1997), i.e., the survey 
region was observed twice with a delay of approximately 
one month between the two observations, and the two 
epochs were then analyzed to find transients. Details 
of the search campaigns can be found in Lidman et al. 
(2005) and Morokuma et al. (2010). 

The Spring 2001 data were processed to find transient 
objects and the most promising candidates were given an 
internal SCP name and a priority. The priority was based 
on a number of factors: the significance of the detection, 
the relative increase in the brightness, the distance from 
the center of the apparent host, the brightness of the 
candidate and the quality of the subtraction. Note how- 
ever that these factors were not applied independently 
of each other, but the priorities were rather based on a 
combination of factors. For example, candidates on core 
were only avoided if they had a small relative brightness 
increase over the span of 1 month. The AGN structure 
function shows that AGNs rarely have strong changes 
over 1 month. 

The candidates discovered in the Spring 2001 campaign 
were distributed to teams working at Keck and Paranal 
observatories for spectroscopic confirmation. The dis- 
tribution was based on the likely redshifts. Candidates 
that were likely to be SNe la at z > 0.7 were sent to 
Keck, while candidates that were thought to be nearer 
were sent to the VLT. In later years, when F0RS2 was 
upgraded with a CCD with increased red sensitivity, the 
most distant candidates would also be sent to the VLT. 
All Subaru 2001 candidates were sent to FOCAS (Faint 
Object Camera And Spectrograph) on Subaru for spec- 
troscopic confirmation (Morokuma et al. 2010). 

In total, four instruments (FORSl on the ESQ VLT, 
ESI and LRIS on Keck, and FOCAS on Subaru) were 
used to determine redshifts and to spectroscopically con- 
firm the SN type. The dates of the spectroscopic runs 
are listed in Table 1 and the observations of individual 
candidates are listed in Table 2 

Only those candidates that were confirmed as SNe la 
were then scheduled for follow-up observations from the 
ground and with HST. In total, six SNe were sent for 
HST follow-up: one from the Subaru 2001 campaign and 
five from the Spring 2001 campaign. The SNe are listed 
in Tables 2 and 3. Finding charts are provided in Fig- 
ure 1. 

Here, we describe the analysis of data that were taken 
with ESI and LRIS. The analyses of the spectra taken 
with FORSl and FOCAS (SN 2001cw, SN 2001go and 
SN 2001gy) are described in Lidman et al. (2005) and 
Morokuma et al. (2010) respectively. The spectra of 
these SN are shown in these papers and will not be re- 
peated here. 

2.1. ESI 

The two highest redshift candidates, 200 Ign and 
2001hb, were observed with the echelette mode of ESI 
(Sheinis et al. 2002). A spectrum taken with the echel- 
lette mode of ESI and the 20" slit covers the 0.39 /um 
to 1.09 /xm wavelength range and is spread over 10 or- 
ders ranging in dispersion from 0.16 A per pixel in the 
bluest order (order 15) to 0.30 A per pixel in the red- 
dest (order 6). The detector is a MIT-Lincoln Labs 



2048 x 4096 CCD with 15 ^im. pixels. The slit width 
was set according to the seeing conditions and varied 
from 0.7" to 1.0", which corresponds to a spectral reso- 
lution of i? 5000. Compared to spectra obtained with 
low-resolution spectrographs, such as FORSl, F0RS2, 
LRIS and FOCAs, the fraction of the ESI spectrum that 
is free of bright night sky lines from the Earth's atmo- 
sphere is much greater. This allows one to de-weight the 
low signal-to-noise regions that overlap these bright lines 
when binning the spectra, a method that becomes in- 
efficient with low-resolution spectra as too much of the 
spectra are de-weighted. Another advantage of ESI was 
that the MIT-LL CCD offered high quantum efficiency at 
red wavelengths and significantly reduced fringing com- 
pared to conventional backside-illuminated CCDs. 

The data were reduced in a standard manner. The bias 
was removed by subtracting the median of the pixel val- 
ues in the overscan regions, the relative gains of the two 
amplifiers were normalized by multiplying one of the out- 
puts with a constant and the data were flat fielded with 
internal lamps. When extracting SN spectra, a bright 
star was used to define the trace along each order, and 
the spectrum of the SN was used to define the center of 
the aperture. Once extracted, the 10 orders were wave- 
length calibrated (using internal arc lamps and cross- 
checking the result with bright OH lines), flux calibrated 
and stitched together to form a continuous spectrum. 

To reduce the impact that residuals from bright OH 
lines have on determining the redshift and classifying the 
candidate, the spectrum was weighted according to the 
inverse square of the error spectrum and then rcbinned 
by a factor of 20, from 0.19 A/pixel to 3.8 A/pixel. The 
binning was chosen so that features from the host were 
not lost. 

The reduced spectra of 2001gn and 2001hb are pre- 
sented in Figures 16 and 17 respectively. 

2.2. LRIS 

The LRIS (Oke et al. 1995) data were taken with the 
400/8500 grating and the GG495 order sorting filter and 
were reduced in a standard manner. The bias was re- 
moved with a bias frame, the pixel-to-pixel variations 
were normalized with flats that were taken with internal 
lamps and the background was subtracted by fltting a low 
order polynomial along detector columns. The fringes 
that were not removed by the flat were removed with a 
fringe map, which was the median of the sky subtracted 
data that was then smoothed with a 5x5 pixel box. The 
spectra were then combined, extracted, and calibrated in 
wavelength and flux. 

The reduced spectrum of 2001gq is shown in Figure 18. 

2.3. Spectral fitting and supernova typing 

Light from the host and the SN are often strongly 
blended in the spectra of high redshift SN. To separate 
the two, we followed the spectral fitting technique de- 
scribed in Howell et al. (2005). 

To classify the SNe, we used the classification scheme 
described in Lidman et al. (2005) and added to it the 
confidence index (CI) described in Howell et al. (2005). 
In the Lidman et al. (2005) scheme, an object is classified 
as a SN la if the Si II features at 4000 A and/or 6150 A 
or the S II W feature can be clearly identified in the 
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Figure 1. SN finding charts. North is up and East is to the left. The SNe are marked with red cross hairs. The images have been created 
by stacking all /-band data taken with the search instrument and all F814W data obtained with WFPC2. The patch widths are l' and 
0.5' for the ground-based and HST images respectively. 
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Table 1 

Lisl of iiislruuK^Lits and 1 (>l(-sc()])t*s used Ibr s|hh;1 i-()sc(>])y. 



Search 



Instr./Tel. 



Detector 



Resolution Observing dates 



Spring 2001 
Spring 2001 
Spring 2001 
Subaru 2001 



FORSl/Antu 
LRIS/Keck I 
ESI/Keck II 
FOCAS/Subaru 



Tektronix 2kx2k CCD 500 

Tektronix 2kx2k CCD 850 

MIT-LL 2kx4k CCD 5000 

SITe 2kx4k CCD 1000 



2001 April 21 - 22 
2001 April 20 
2001 April 21 - 24 
2001 May 26 - 27 



Table 2 

Summary of the spectroscopic observations. For completeness, the observations that are reported in 
Lidman et al. (2005) and Morokuma et al. (2010) are also included. The Galactic extinction, E{B — V), 
from Schlegel et al. (1998), is presented together with the coordinates for each SN. 







a (2000) 




5(2000) 


E{B - V) 


MJD 


Instrument and 


Exp. 


lAU name 


Search 










(days) 


telescope 


(s) 


SN 2001CW 


Subaru 01 


15''23™06^ 


.3 


+29°39'32" 


0.024 


52056.6 


FOCAS/Subaru* 


4200 


SN 2001gn 


Spring 01 


14h01™59= 


.9 


+05°05'00" 


0.028 


52023.1 


ESI/Kcck II'' 


9700 


SN 2001go 


Spring 01 


141^02^00= 


.9 


+05°00'59" 


0.027 


52021.3 


FORSl/Antu'' 


2400 


SN 2001gq 


Spring 01 


14h01™51= 


.4 


+04°53'12" 


0.027 


52020.3 


LRIS/Kcck^' 


3600 


SN 200 Igy 


Spring 01 


13'^57™04= 


.5 


+04°31'00" 


0.030 


52021.3 


FORSl/Antu'' 


2400 


SN 2001hb 


Spring 01 


13''57™11= 


.9 


+04°20'27" 


0.032 


52024.3 


ESI/Keck II'' 


3600 



Long slit 
Echellette 
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spectrum or if the spectrum is best fit with the spectra 
of nearby SN la and other types do not provide a good 
fit. We quahfy the classification with the keys "Si II" or 
"SF" in column 3 of Table 3 depending on whether the 
classification was done by identifying features or by using 
the fit. In the scheme described in Howell et al. (2005), 
these SNc would have a CI of 5 and 4, respectively. 

Less secure candidates are classified as la*. The as- 
terisk indicates some degree of uncertainty. Usually, this 
means that we can find an acceptable match with nearby 
SNe la; however, other types, such as SNe Ibc, also re- 
sult in acceptable matches. These SNe have a confidence 
index of 3. 

Rcclshifts based on the host have an accuracy that is 
better than 0.001, and are, therefore, quoted to three 
decimal places. Redshifts based on the fit are less accu- 
rate. For completeness, the redshifts and classifications 
reported in Lidman et al. (2005) and Morokuma et al. 
(2010) are also included. 

The agreement between the phase, tspoc, of the best 
fit template and the corresponding phase, thc, obtained 
from the light curve fit is also shown in Table 3. The 
weighted average difference for all six spectra is At = 
—0.4 days with a dispersion of 2.0 days. The dispersion 
is similar in magnitude to that found in other surveys 
(Hook et al. 2005; Foley et al. 2008). 

3. PHOTOMETRIC OBSERVATIONS 

A total of nine different instruments, listed in Table 4, 
were used for the photometric follow-up of the SNc de- 
scribed in this work. 

All observations arc listed in Tabic 12. Here the Modi- 
fied Julian Date (MJD) is the weighted average of all im- 
ages taken during a given night except for the NIR data 
where data taken over several nights were combined. We 
do not report the MJD for combined reference images 
that were taken over several months. 

3.1. Ground-based optical observations and reductions 

We obtained ground-based optical follow-up data of 
the SNc through different combinations of passbands, 
shown in Figure 2, similar to Besscl R and / (Bessell 
1990), and SDSS i (Fukugita et al. 1996). Here we 
used FORSl (Appenzeller et al. 1998) at VLT and SuSI2 
(D'Odorico et al. 1998) at NTT in addition to the 
search instruments. The SuSI2, MOSAIC II, CFH12k 
and SuprimeCam data were obtained in visitor mode, 
while the FORSl observations were carried out in ser- 
vice mode. 

All optical ground data were reduced (Raux 2003) in a 

standard manner including bias siibtraction, flat fielding 
and fringe map subtraction using the IRAF^'^ software. 

3.2. Ground-based IR observations and reduction 

The two most distant SNe in the sample, 2001hb and 
2001gn, were also observed from the ground in the near- 
IR. Both NIRI (Hodapp et al. 2003) and ISAAC (Moor- 
wood et al. 1999) were used to observe 2001hb, while 
2001gn was observed with ISAAC only. 

^® IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under the cooperative agreement 
with the National Science Foundation. 



The ISAAC observations were done with the Js filter 
and the NIRI observations were carried out with the J 
filter. The transmission curves of the filters are similar 
to each other, and the transmission curve of the latter 
is shown in Figure 2. The red edges of the filters are 
defined by the filters and not by the broad telluric ab- 
sorption band that lies between the J and H windows, 
and the central wavelength is slightly redder than the 
central wavelength of the J filter of Persson et al. (1998). 
Compared to traditional J band filters, photometry with 
the ISAAC Js and NIRI J band filters is less affected by 
water vapor and is therefore more stable. 

The ISAAC observations were done in service mode 
and the data were taken on 14 separate nights, starting 
on 2001 May 7 and ending on 2003 May 30. Individual 
exposures lasted 30 to 40 s, and three to four of these 
were averaged to form a single image. Between images, 
the telescope was offset by 10" to 30" in a semi-random 
manner, and typically 20 to 25 images were taken in this 
way in a single observing block. The observing block was 
repeated several times until sufficient depth was reached. 

The data, including the calibrations, were first pro- 
cessed to remove two electronic artifacts. In about 10% 
of the data, a difference in the relative level of odd and 
even columns can be seen. The relative difference is a 
function of the average count level and it evolves with 
time, so it cannot be removed with flat flelds. In those 
cases where the effect is present, the data axe processed 
with the eclipse^^ odd-even routine. The second arti- 
fact, an electronic ghost, which is most easily seen when 
there are bright stars in the field of view, is removed with 
the eclipse ghost routine. 

The NIRI observations were done in queue mode and 
the data were taken on 4 separate nights, starting on 2001 
May 25 and ending 2002 Aug 5. Individual exposures 
lasted 60 s. Between images, the telescope was offset 
by 10" to 30" in a semi-random manner, and typically 
60 images were taken in this way in a single observing 
sequence. The sequence was repeated several times until 
sufficient depth was reached. 

Both the ISAAC and NIRI data were then reduced in 
a standard way with the IRAF XDIMSUM package and 
our own IRAF scripts. From each image, the zero-level 
offset was removed, a fiatfield correction was applied, and 
an estimate of the sky from other images in the sequence 
was subtracted. Images were then combined with in- 
dividual weights that depend on the median sky back- 
ground and the image quality. 

3.3. HST observations and reduction 

Observing SN la at high-2 from space has an enormous 
advantage for accurately following their light curves. The 

absence of the atmosphere and the high spatial resolution 
allows high signal-to-noise measurements. The high spa- 
tial resolution also helps minimize host contamination 
through focusing the light over a smaller area. Space 
also permits observations at longer wavelengths where 
the limited atmospheric transmission and the high back- 
ground degrade ground-based data. 

High quality follow-up data were obtained for all six 
SNe using the Wide Field Planetary Camera 2 (WFPC2) 
on the Hubble Space Telescope (HST) during Cycle 9. 

http: / / www.eso.org/ projects / aot / eclipse / 
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Table 3 

Classifications, redshifts, confidence index (Howell et al. 2005) and matches to nearby SN 
spectra (or templates). Unless indicated otherwise, the uncertainty in the redshift is 

0.001. Here tspoc shows the spectroscopic phase relative to maximum light in the B-band 
of the template, (lc shows the rest frame phase of the observation relative to the fitted 
B-band maximum (MJDmax in Table 6), and At is the difference. The confidence index 
(CI) and spectroscopic classification are explained in detail in the text. We also specify 
whether the classification was done from the fit ("SF") or if "Si H" could be identified. 



lAU name 


CI 


Spectroscopic 

classification 


Redshift 


Template 

match 


^Spec 

(days) 


(days) 


At 

(days) 


2001CW 


3 


la* 


SF 


0.95 ± 0.01 


1989B 


-5 


-5.1 ±0.4 


+0.1 


2001gn 


3 


la* 


SF 


1.124 


1990N 


-7 


-3.4 ±0.1 


-3.6 


2001go 


5 


la 


Si 11 


0.552 


1992A 


+5 


+6.5 ±0.2 


-1.5 


2001gq 


3 


la* 


SF 


0.672 


1999bp 


-2 


-5.1 ±0.2 


+3.1 


2001gy 


5 


la 


Si n 


0.511 


1990N 


-7 


-6.9 ±0.2 


-0.1 


2001hb 


3 


la* 


SF 


1.03 ± 0.02 


1989B 


-7 


-6.8 ±0.5 


-0.2 



Table 4 

Instruments and telescopes used for photometric follow-up. For WFPC2 we only list the 
properties of the PC camera which was always used for targeting the SNe. 



Tel./Instr. 


Scale 

C'/Px) 


FOV 

(') 


Detectors 


Search 


CFHT/CFH12k 


0.206 


42x28 


12 MIT 2kx4k CCD 


Spring 2001 


CTIO/MOSAIC H 


0.27 


36x36 


8 SITe 2kx2k CCD 


Spring 2001 


VLT/FORSl 


0.20 


6.8x6.8 


1 Tektronix 2kx2k CCD 


Spring 2001 


VLT/ISAAC 


0.1484 


2.5x2.5 


Hawaii Ikxlk HgCdTe 

array 


Spring 2001 


NTT/SuSI2 


0.08 


5.5x5.5 


2 EEV 2kx4k CCD 


Spring 2001 


Gemini/NIRI 


0.1171 


2.0x2.0 


Aladdin Ikxlk InSb array 


Spring 2001 


HST/ACS 


0.05 


2.4x2.4 


2 SITe 2kx4k CCD 


Spring 2001 


HST/WFPC2 (PC) 


0.046 


0.61x0.61 


1 Loral Aerospace 

800x800 CCD 


Both 


Subaru / SuprimeCam 


0.20 


34x27 


10 MIT/LL 2kx4k CCD 


Subaru 2001 



All but two (2001cw and 2001gq) of the objects also have 
SN-free reference images taken during Cycle 10, with the 
Advanced Camera for Survey (ACS). WFPC2 consists 
of four 800x800 pixel chips of which one, the Planetary 
Camera (PC), has twice the resolution of the others. The 
SNe were always targeted with the lower left corner of the 
PC in order to place them closer to the readout amplifier 
so that the effects of charge-transfer inefficiency would be 
reduced. Images of the same target were obtained with 
roughly the same rotation for all epochs. 

Each SN was followed in two bands. All SNe were ob- 
served in the F814W filter. Additionally, the F675W fil- 
ter was used for the three SNe at 2; < 0.7 and the F850LP 
filter was used for the high redshift targets. These filters 
were chosen to match the filters used on the ground and 
correspond approximately to rest frame UBV hands as 
illustrated in Figure 2. 

The data were reduced with software provided by the 



Space Telescope Science Institute (STScI). WFPC2 im- 
ages were processed through the STScI pipeline and then 

combined for each epoch to reject cosmic rays using the 
crrej task which is part of ST SDAS'^" IRAF package. 

The ACS images were processed using the multi-drizzle 
(Fruchter & Hook 2002) software, which also corrects for 
the severe geometric distortion of the instrument. For 
this we used the updated distortion coefficients from the 
ACS Data Handbook in November 2006, and drizzled the 
images to the resolution of the WFPC2 images 0.046". 
Note that the WFPC2 images were not corrected for ge- 
ometric distortion at this stage. 

4. PHOTOMETRY OF THE GROUND-BASED DATA 

The Space Telescope Science Data Analysis System (STS- 
DAS) is a software package for reducing and analyzing astronom- 
ical data. It provides general-purpose tools for astronomical data 
analysis as well as routines specifically designed for HST data. 
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Figure 2. Illustration of filters (colored) used in the SN campaign together with the rest frame Bossel filters (gray). The red (ground- 
based) and blue (HST) filters were used for the observations of the three SNe below (upper panel) and above (lower panel) z = 0.7 
respectively. The filters have been blueshifted to the mean of the two sub samples. Note that the Bessel filters shown here differ slightly 
from the corresponding filters used at the different telescopes. Not shown in the figure is the SuprimeCam i filter, which is close to the I 
filter and was only used for 2001cw. We show the J band filter curve for NIRl which is very similar to the ISAAC Js filter. The atmospheric 
transmission for 2 mm of precipitable water vapor, typical for Paranal and Mauna Kea, is also plotted between 10000 — I6OOOA in the 
observer frame. The transmission curve was provided by Alain Smette (private communication). 



The photometry technique apphed to the optical 
ground-based data is the same as the one apphed in 
Amanuhah et aL (2008), except for the SuprimeCam i 
band data. This is also very similar to the method (Fab- 
bro 2001) used in Astier et al. (2006), and is briefly sum- 
marized here: 

1. Each exposure of a given SN in a given passband 
was aligned to the best seeing photometric refer- 
ence image. 



In order to properly compare images of different 
image quality we fitted convolution kernels, Ki, 
modelled by a linear decomposition of Gaussian 
and polynomial basis functions (Alard & Lupton 
1998; Alard 2000), between the photometric refer- 
ence and each of the remaining images, i, for the 
given passband. The kernels were fitted by using 
image patches centered on fiducial objects across 
the field. 



The background sky level for each image, i, is not 
expected to have any spatial variation for a small 
patch, Ii{x,y), centered on the SN with a radius 
of the worst seeing FWHM. We assume that the 
patch can be modelled by a point spread function 
(PSF) at the location of the SN, a model of the 



host galaxy and a constant offset for the sky, 

h{x, y) = h- [K, (g) PSF] (x ~XQ,y- yo) + 
+ [K, (g) G] (x, y) +S^. 

For a time-series of such patches we simultaneously 
fit the SN position, (a;o,j/o), and brightnesses, /i, 
host model, G{x, y), and background sky levels. Si. 
We use a non-analytic host model with one param- 
eter per pixel. This means that the model will be 
degenerate with the SN and the sky background. 
We break these degeneracies by fixing the SN flux 
to zero for all reference images and fix the sky level 
to zero in one of the images. Figure 3 shows an ex- 
ample of image patches, galaxy model residuals and 
resulting residuals when the full model has been 
subtracted, for increasing epochs. 

The SN light curves were obtained in this manner for one 
filter and instrument at a time where SN-free reference 
images were available. In the cases where SN-free refer- 
ences did not exist for a given telescope, reference images 
obtained with other telescopes were used instead. 

For the IR image of 2001hb, we performed aperture 
photometry directly on the images, assuming the SN to 
be hostless for the purpose of J band photometry, since 
no host could be detected at the limit of the ACS refer- 
ences (see below). For the IR image of 2001gn, we carried 
out steps 1-2 above, and then subtracted the reference 
image (Figure 4) and did aperture photometry on the 
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Figure 3. Ground-based /-band images (patches) of 2001gq (three left columns; widtli 6.0") and 2001go (three right columns; widtii 
4.7"). Each triplet represents from left to right tiic fully reduced data, the galaxy model subtracted data, and a profile plot where the 
full galaxy -|- PSF model has been subtracted from the data. The profile plot shows the deviation from zero in standard deviations of the 
noise versus the distance from the SN position in pixels. Epoch increases from top to bottom. The images in the last row are the reference 
image, obtained approximately one year after discovery. The different image sizes (patches) reflect differences in image quality. A small 
image represents good seeing conditions. 
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resulting imago. In both cases the SN positions from the 
HST images were used for ccntroiding the aperture and 
the diameter was chosen to maximize the signal-to-noise 
ratio. The fluxes are corrected to larger apertures by 
analyzing bright stars in the same image. 

The i band data for 2001cw was analyzed together with 
a larger sample of SNc discovered at Subaru. The details 
of this analysis are given in Yasuda et al. (in prepara- 
tion). 

4.1. Calibration of the ground-based data 

Nightly observations of standard stars wore not avail- 
able for all optical instruments and filters. Instead 
the recipe from AmanuUah et al. (2008), using SDSS 
(Adelman-McCarthy et al. 2007) measurements of the 
field stars was applied. However, the SDSS filter sys- 
tem (Fukugita et al. 1996) differs significantly from the 
filters used in this search except for the SuprimeCam 
i band. To overcome this difference, we fitted rela- 
tions between the SDSS filter system and the Landolt 
system (Landolt 1992) in a similar manner to Lupton 
(2005), using stars with SDSS and Stetson photometry 
(Stetson 2000). Stetson has been publishing photome- 
try of a growing list of faint stars that is tied to the 
Landolt system within 0.01 mag (Stetson 2000). The 
most up to date version can be obtained from the Cana- 
dian Astronomy Data Center^^. In contrast to Lupton 
(2005), we applied magnitude, uncertainty and color cuts 
(15 < rsDSS < 18.5, am < 0.03 mag in tsdss and isDSS 
and —0.5 < tsdss — isuss < 1-0) for the stars that went 
into the fit. We also applied a 5a outlier cut after our 
initial fit (and lost about ~ 6% of the sample). Af- 
ter refitting, the following relations were derived for the 
Landolt R and / filters 

R - rsDSS = -0.13 - 0.32 • (rsDSS - «SDSs) 

I - isDss = -0.38 - 0.26 ■ (rsDss - «sdss), 

which are also shown in Figure 5. The results for R — 
''SDSS are close to the ones derived by Lupton (2005) as 
well as to Tonry et al. (2003) who performed a similar 
operation. Neither Tonry et al. (2003) nor Lupton (2005) 
present fits for / - isdss vs rsoss - «Sdss- 

By forcing x^/dof = 1, we determined that there is a 
scatter of 0.03 mag coming from intrinsic spectral distri- 
butions. This is a systematic uncertainty for individual 
stars, but will average out when a big sample is used, as- 
suming that the color distribution of the sample is similar 
to the stars used to derive the relations. 

The transformations were applied to the SDSS stellar 
photometry of our SN field, and these were then used as 
tertiary standard stars in order to tic the SN photometry 
to the Landolt system. The flux of the stars was deter- 
mined on the photometric reference for each light curve 
build using the same method as for the SN and with the 
same PSF model that was used for fitting the SN fluxes. 
A zero-point relation of the form 

m + 2.51ogio/ = ZP + cx-(ii-/) (1) 

could then be fitted between the measured stellar fluxes, 
/, and their Landolt magnitudes, m. Here ZP is the zero 
point and cx is the color term for the filter. We applied 

http : //mra . cade . hia . nr c . gc . ca/community/STETSON/ archive/ 



the same color cuts (—0.5 < rgDss — *sdss < 1-0) to the 
stars that went into the fit. 

Unfortunately we did not have enough stars over a wide 
enough color range to accurately fit the color term. In- 
stead, we used values from the literature or from the ob- 
servatories, which are summarized in Table 5. The zero- 
points could then be derived from equation (1). The val- 
ues obtained this way are the sums of three components; 
the instrumental zero-point, the aperture correction to 
the PSF normalization radius and the atmospheric ex- 
tinction for the given airmass, and are shown along with 
the SN fluxes in Table 12. 

We also calculated color terms synthetically by using 
Landolt standard stars that have extensive spectropho- 
tometry from Stritzinger et al. (2005). The synthetic 
Vega magnitudes of the stars were calculated by multi- 
plying the spectra and the Vega spectrum (Bohlin 2007) 
with the filter and instrument throughputs provided by 
the different observatories. The color terms were then 
fitted by assuming a linear relation for the deviation be- 
tween the synthetic and the Landolt magnitudes as a 
function of the Landolt color. The resulting fitted syn- 
thetic terms are also presented in Table 5. 

The difference is < 0.03 mag for all but the CFH12k 
/-band, where there is a significant discrepancy. A mis- 
match between the effective filter transmission curves we 
use and the true photometric system is most likely the 
origin of this. We have investigated if the deviation could 
be explained by the differences in quantum efficiencies 
between the different chips of the detector which are pre- 
sented on the CFH12k webpage. However these differ- 
ences propagated to the color term results in relatively 
low scatter and deviates significantly from the measured 
value. 

For the purpose of fitting the zero-points, we can use 
the measured color terms, but for light curve fitting er- 
roneous filter transmissions could introduce systematic 
effects. In order to study the potential impact on the 
light curve fits we modified the filters to match the mea- 
sured color terms. The modifications were implemented 
by either shifting or clipping the filters until the synthetic 
color terms matched all measured values for a given fil- 
ter (e.g. R OY I vs R — I, V — I and _B — / for the 
ELIXIR measurements for the /-band) . The light curves 
were then refltted using the modifled filter transmissions 
and the resulting values were compared. Since the light 
curves are tightly constrained by the high precision HST 
data, the modified filter only leads to negligible differ- 
ences (less than 30 % of the statistical uncertainty) in 
the fitted parameters. 

The J-band data, on the other hand, were calibrated 
using G-type standard stars from the Persson LCO stan- 
dard star catalogue (Persson et al. 1998). Since the 
ISAAC and NIRI J-band filters are slightly redder and 
narrower than the Persson J-band filter, we subtracted 
0.012 magnitudes from the ISAAC and NIRI zero points 
to place the ISAAC and NIRI IR photometry onto the 
natural system. Since the data were taken over many 
nights, the photometry was carefully cross checked. Dif- 
ferences in the absolute photometry visually amount to 
less than 0.02 magnitudes, which we conservatively adopt 
as our zero-point uncertainty. 

5. PHOTOMETRY OF THE HST DATA 
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Figure 4. VLT/ISAAC J band observations of the distant SN la, 2001gn (North is up and East is to the left). On the left, we show a 27" 
wide and 51" high image with the SN; in the middle we show the reference image, which was taJien 2 years later. The images on the right 
is the subtraction between the two and shows the SN with a signal-to-noise ratio of 10. Each image is the result of 10 hours of integration 
in good conditions. The image quality in the left hand image is FWHM=0.39". 







Am 








Am 



Figure 5. Stars observed both by Peter Stetson (Stetson 2000) and SDSS for which —0.5 < rgDSS — *SDSS < 1-0 together with fitted 
linear relations between the SDSS rgDSS, *SDSS and the Landolt R, I magnitudes. The left column shows the relation as a function of color, 
while the middle column shows the residuals once the fitted relation has been subtracted, and a histogram of the residuals are presented 
in the right column. 



A modified version of tlie pliotonictric teclmiquc from 
Knop et al. (2003) was used for the HST data presented 
here. This is similar to the method used for the ground- 
based optical data above, but instead of aligning and 
resampling all images to a common frame, the host+SN 
model is resampled to each individual image. A proce- 
dure like this is preferred when the PSF FWHM is of the 
same order as the pixel scale, and it also preserves the 
image noise properties. 



Linear geometric transformations were first fitted from 
each image, fc, in a given filter to the deepest image of the 
field using field objects. Due to the similar orientation 
of the WFPC2 images, linear transformations were suffi- 
cient for these, and the geometric distortion of WFPC2 
could be ignored. This was however not the case for 
the transformations between the ACS and WFPC2 im- 
ages and the distortion was then hardcoded into the fit- 
ting procedure. Due to the sparse number of objects 
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Table 5 

Measured and synthetic color terms for the three 
instruments that were used as photometric references. 





Measured 


Synthetic 


Inst./Tel. 




C/ 






CFH12k/CFHT 


0.031=- 


0.107=- 


0.066 


-0.023 


MOSAIC II/CTIO 




O.OSO'' 




0.007 


FORSl/VLT 


0.034"= 


-0.050"= 


0.040 


-0.076 



=■ http: //www. cfht .hawaii .edu/Instruments/Elixir/f liters .html 

Miknaitis et al. (2007) 

"= http://www.eso.org/observing/dfo/quality/FORSl/qc/photcoeff/photcoeffs_forsl.html 



in the tiny PC field the accuracy of all transformations 
were only good to < 1 pixel (~ 0.5 FWHM). Unfortu- 
nately, using objects from the remaining three chips for 
the alignment did not lead to improved accuracy, which 
is probably due to small movements of the chips between 
exposures (Anderson & King 2003). This alignment pre- 
cision was not enough for PSF photometry, and we there- 
fore allowed the SN position, (.t/jo, y/co) to float for the 
individual images which increased the alignment preci- 
sion by a factor of 10. Allowing this extra degree of 
freedom could bias the results toward higher fluxes since 
the fit will favor positive noise fluctuations. However, as 
in Knop et al. (2003), this was shown to be of minor im- 
portance by studying the covariance between the fitted 
flux and SN position. 

The full model used to describe each image patch can 
be expressed as 

h{xk,yk) = fk ■ PSFk{xk - Xko,yk - yko)+ , , 
+ G{xk,yk,aj) + Sk ■ 

Here Ik is the value in pixel (xk,yk) on image k, fk is 
the SN flux, PSF;; the point spread function, G the host 
galaxy model that is parameterized by aj and Sk the 
local sky background. The fits were carried out using 
a minimization approach using MINUIT (James & 
Roos 1975). 

Four SNc in the sample had ACS reference images, and 
for these cases we used field objects to fit non-analytic 
convolution kernels, K, between the PC chip and the 
drizzled ACS image. These accounted for the difference 
in quantum efficiency and PSF shape between the two 
instruments. When kernels were used, the uncertainties 
of individual pixels were propagated and the correlation 
between pixels produced by the convolution were taken 
into account. Also, in this case equation (2) above was 
modified so with both the PC images and the PC PSF 
being convolved with the fitted kernel. 

The PSF, PSFfe, of the WFPC2 PC chip was simu- 
lated for each filter and pixel position using the Tiny Tim 
software (Krist 2001) and normalized to the WFPC2 cal- 
ibration radius 0.5". We also did an extensive test where 
we iteratively updated the simulated PSF based on the 
knowledge of the SN epoch, and therefore the approxi- 
mate spectral energy distribution (SED), but this did not 
have any significant effect on the fitted fluxes. The Tiny 
Tim PSF was generated to be subsampled by a factor of 
10. For each iteration in the fitting procedure, any shift 
of the PSF position was first applied in the subsampled 



space. The PSF was then re-binned to normal sampling 
and convolved with a charge diffusion kernel (Krist 2001) 
before it was added to the patch model. 

For three of the SNe, 2001cw, 2001gq and 2001hb 
we used an analytical model for the host galaxy. Both 
2001cw and 2001gq are offset from the core of their re- 
spective host galaxies and we therefore chose not to ob- 
tain SN-free references images for these. Instead the 
hosts could be modeled by a second order polynomial 
and constrained by the galaxy light in the vicinity of the 
SNe. For 2001hb, we did obtain a deep ACS reference 
image but no host could be detected. For this SN we 
used a simple plane to model the host which was further 
constrained by fixing the SN flux to zero for the ACS 
reference. One caveat with analytical host modelling in 
general is that the patch size must be chosen with care. 
The model will only work if there is no dramatic change 
in the background across the patch, which is an assump- 
tion that is likely to fail if the patch is too large and in- 
cludes the host galaxy core. On the other hand the patch 
can not be too small either in order to successfully break 
the degeneracy between the SN and the background. 

To make sure that the choice of host model does not 
bias the fitted SN fluxes we required, in addition to clean 
residuals once the SN + host was subtracted from the 
data, that the fitted SN fluxes were insensitive to vari- 
ations in the patch size of a few pixels. Further, for 
2001gq, we tested the host model by putting fake SNe 
at different positions around the core of the host. The 
distances to the core and the fluxes of the fakes were al- 
ways chosen to match the corresponding values for the 
real SN, and the retrieved photometry was always within 
the expected statistical uncertainty. 

The recipe described above could not be used for the 
remaining three SNe in the sample, since they were lo- 
cated too close to the cores of their host galaxies. Addi- 
tionally, the hosts have small angular sizes and the light 
gradients in the vicinity of the SNe were steep enough to 
lead to biased photometry due to the coarse geometric 
alignment. To overcome this we changed the photometric 
procedure slightly. Instead of fitting the SN position on 
each image, we chose to fit it only on the geometric refer- 
ence image, and then introduce a free shift for the whole 
model. That is, in this case we used the galaxy model -|- 
SN for the patch alignment. However, this procedure did 
force us to apply some constraints on the galaxy mod- 
elling. Using a non-analytic pixel model, the approach 
used for the ground-based optical data, was not feasi- 
ble. It slowed down the fits considerably and rarely con- 
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verged. Instead we chose to use the ACS images of the 
host galaxies directly. The ACS references were much 
deeper then the WFPC2 data and choosing this proce- 
dure did not increase the uncertainty of the fitted fluxes. 

A general problem with doing photometry on HST im- 
ages is that CCD photometry of faint objects over a low 
background suffers from an imperfect charge transfer, 
which will lead to an underestimate of the flux. We 
used the Charge Transfer Efficiency (CTE) recipe for 
point sources from Dolphin (2009). The correction for 
our data is usually around 5 %-8 % but it can be as large 
as ^ 17%. The uncertainties of the corrections were 
propagated to the flux uncertainties. The corrected SN 
fluxes arc given in Table 12 together with the instrumen- 
tal zero-points, which were also obtained from Andrew 
Dolphin's webpage 

6. LIGHT CURVE FITTING 

SN la that have bluer colors or broader light curves 
tend to be intrinsically brighter (Phillips 1993; Tripp 
1998). Several methods of combining this information 
into an accurate measure of the relative distance have 
been used (Riess et al. 1996; Goldhaber et al. 2001; Wang 
et al. 2003; Guy et al. 2005, 2007; Jha et al. 2007; Conley 
et al. 2008). 

K08 consistently fitted all light curves using the SALT 
(Guy et al. 2005) fitter, which is built on the SN la SED 
from Nugent et al. (2002). In this paper, we use SALT2 
(Guy et al. 2007), which is based on more data. 

Conley et al. (2008) compared the performances of dif- 
ferent light curve fitters while also introducing their own 
empirical fitter, SiFTO, and concluded that SALT2 along 
with SiFTO perform better than both SALT (which is 
conceptually different from its successor SALT2) and 
MLCS2k2 (Jha et al. 2007) when judged by the scat- 
ter aroimd the best-fit luminosity distance relationship. 
Furthermore, SALT2 and SiFTO produce consistent cos- 
mological results when both are trained on the same 
data. Recently KS09 made a thorough comparison be- 
tween SALT2 and their modified version of MLCS2k2 
(Jha et al. 2007) for a compilation of public data sets, in- 
cluding the one from the SDSS SN survey. The two light 
curve fitters result in an estimate of w (for a flat wCDM 
cosmology) that differs by 0.2. The difference exceeds 
their statistical and systematic (from other sources) er- 
ror budgets. They determine that this deviation origi- 
nates almost exclusively from the difference between the 
two fitters in the rest-frame J7-band region, and the color 
prior used in MLCS2k2. They also noted that MLCS2k2 
is less accurate at predicting the rest- frame U -band using 
data from filters at longer wavelengths. 

This difference in J7-band performance is not surpris- 
ing: observations carried out in the observer-frame U- 
band arc in general associated with a high level of uncer- 
tainty due to atmospheric variations. While the train- 
ing of MLCS2k2 is exclusively based on observations of 
nearby SNe, the SiFTO and SALT2 training address this 
difficulty by also including high redshift data where the 
rest-frame J7-band is observed at redder wavelengths. 
This approach also allows these fitters to extend blue- 
ward of the rest-frame J7-band. 

In addition, for this paper, we have conducted our own 
test validating the performance of SALT2 by carrying out 
the Monte-Carlo simulation described in 7.3.6, where we 



compare the fitted SALT2 parameters to the correspond- 
ing real values for mock samples with poor cadence and 
low signal-to-noise drawn from individual well-measured 
nearby SNe. 

Given these tests that have been carried out on SALT2, 
and its high redshift source for rest-frame ?7-band, we 
have chosen to use SALT2 in this paper. 

6.1. SALT2 

The SALT2 SED model has been derived through a 
pseudo-Principal component analysis based on both pho- 
tometric and spectroscopic data. Most of these data 
come from nearby SN la data, but SNLS supernovae 
are also included. To summarize, the SALT2 SED, 
F(SN,p, A), is a function of both wavelength. A, and 
time since B-band maximum, p. It consists of three com- 
ponents; a model of the time dependent average SN la 
SED, Afo(p, A), a model of the variation from the aver- 
age, Mi{p, A), and a wavelength dependent function that 
warps the model, CL{X). The three components have 
been determined from the training process (Guy et al. 
2007) and are combined as 

F(SN,p,A)=xo X [Moip,X) + xi x Mi(p,A)] x 

X exp [c x CL{X)] , 

where xq, xi and c are free parameters that are fit for 
each individual SN. 

Here, xq, describes the overall SED normalization, xi, 
the deviation from the average decline rate (xi = 0) 
of a SN la, and c, the deviation from the mean SN la 
B — V color at the time of i3-band maximum. These 
parameters are determined for each observed SN by fit- 
ting the model to the available data. The fit is carried 
out in the observer frame by redshifting the model, cor- 
recting for Milky Way extinction (using the CCM-law 
from Cardelli et al. (1989) with Rv = 3.1), and multi- 
plying by the effective filter transmission functions pro- 
vided by the different observatories. All synthetic pho- 
tometry is carried out in the Vega system using the 
spectrum from Bohlin (2007). Following Astier et al. 
(2006) we adopt the magnitudes {U,B,V,Rc,Ic) = 
(0.020,0.030,0.030,0.030,0.024) mag (Fukugita et al. 
1996) for Vega. For the near-infrared we adopt the values 
J = and iJ = 0. 

In the fit to our data, we take into account the cor- 
relations introduced between different light curve points 
from using the same host galaxy model. We also chose 
to run SALT2 in the mode where tlie diagonal of the 
covariance matrix is updated iteratively in order to take 
model and /^-correction uncertainties into account. See 
Guy et al. (2007) for details on this. The Milky Way red- 
dening for our supernovae from the Schlegel et al. (1998) 
dust maps is given in Table 2. The results of the fits are 
shown in Table 6 and plotted in Figure 6 together with 
the data. 

The three parameters 



^max ^ _2.51ogi 



i^(SN,0,A) AdA 



, xi and c 



can for each SN be combined to form the distance mod- 
ulus (Guy et al. 2007), 

MB = ml"" - Mb = m'S^ + a ■ xi - 13 ■ c - Mb , (3) 
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Figure 6. SALT2 light curve fits togettier witli tlie data for the six SNe presented here. Light curves for different bands/instruments, 
indicated by the color coding, have been offset. The dotted lines following the curves represents the model errors. 
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where Mb is the absolute _B-band magnitude. The re- 
sulting color and light curve shape corrected peak B- 
band magnitudes, rrf£" , are presented in the third col- 
umn of Table 6. The parameters a, /3 and Mb are nui- 
sance parameters which are fitted simultaneously with 
the cosmological parameters. 

7. THE Union2 COMPILATION 

K08 presented an analysis framework for combining 
different SN la data sets in a consistent manner. Since 
then two other groups (H09 and KS09) have made sim- 
ilar compilations, using different fitters. In this work 
we carry out an improved analysis, using and refining 
the approach of K08. We extend the sample with the 
six SNe presented here, the SNe from Amanullah et al. 
(2008), the low- 2; and intermediate- 2; data from Hicken 
et al. (2009a) and Hohzman et al. (2008) respectively''2_ 

First, all light curves are fitted using a single light curve 
fitter (the SALT2 method) in order to eliminate differ- 
ences that arise from using different fitters. For all SNe 
going into the analysis we require: 

1. data from at least two bands with rest-frame cen- 
tral wavelengths between 2900A and 7000A, the 
default wavelength range of SALT2 

2. that there is at least one point between —15 days 
and 6 rest frame days relative to the B-band max- 
imum. 

3. that there are in total at least five valid data points 
available. 

4. that the fitted xi values, including the fitted un- 
certainties, lie between — 5 < xi < 5. This is a 
more conservative cut than that used in K08 and 
results in several poorly measured SNe being ex- 
cluded. Part of the discrepancy observed by KS09 
when using different light curve models could be 
traced to poorly measured SNe. 

5. that the CMB-centric redshift is greater than z > 
0.015. 

We also exclude one SN from the Union compilation that 

is 1991bg-like, which neither the SALT nor the SALT2 
models are trained to handle. Note that another 1991bg- 
like SN from the Union compilation was removed by the 
outlier rejection. All SNe la considered in this compila- 
tion are listed in Table 13. For each SN, the redshift and 
fitted light curve parameters are presented as well as the 
failed cuts, if any. 

ft should be pointed out that the choice of light curve 
model also has an impact on the sample size. Using 
SALT2 will allow more SNe to pass the cuts above, since 
the SALT2 model covers a broader wavelength range 
than SALT. This is particularly important for high-2; 
data that heavily rely on rest-frame UV data. For ex- 
ample, two net SNe would have been cut from the Riess 
et al. (2007) sample with the SALT model. 

The SALT2 fit results for these samples are pre- 
sented along with the entire Union2 compilation fits at 
http : //supernova . Ibl . gov /Union/. 



7.1. Revised HST zero-points and filter curves 

Since Riess et al. (2007), the reported zero-points of 

both NICMOS and ACS were revised. For the FllOW 
and F160W filters of NICMOS, the revision is substan- 
tial. Using the latest calibration (Thatte et al. 2009, and 
references within), the revised zero-points are, for both 
filters, approximately 5 % fainter than those reported in 
Riess et al. (2007) and subsequently used by K08. 

For SNe la at 2: > 1.1, observations with NICMOS 
cover the rest frame optical, so the fitted peak S-band 
magnitudes and colors and the corrected i?-band mag- 
nitudes of these SNe la depend directly on the accuracy 
of the NICMOS photometry. With the new zero-points, 
SNe la at z > 1.1 are measured to be fainter and bluer. 
Our current analysis also corrects an error in the NIC- 
MOS filter curves that were used in K08, which also acts 
in the same direction. 

In the introduction, we had noted that almost all SNe 
at z > 1.1 were redder than the average SN color over the 
redshift interval 0.3 to 1.1. This is surprising as redder 
SNe are also fainter and should therefore be the harder 
to detect in magnitude limited surveys. K08 noted that 
these SNe, after light curve shape and color corrections, 
are also on average ^ 0.1 mag brighter than the line 
tracing the best fit ACDM cosmology. They also noted 
that this was the reason for the relatively high value for 
the binned value of w in the [0.5, 1] redshift bin. 

After taking the NICMOS zero-point and filter updates 
just discussed into account, we repeated the original K08 
analysis. This made the NICMOS observed SNe up to ~ 
0.1 mag fainter, and there no longer is a significant offset 
from the best fit cosmology. Nor are these SNe unusually 
red when compared to SNe over the redshift interval 0.3 
to 1.1. For SALT2 the SNe at z > 1.1 have an average 
color of c = 0.06 ± 0.03, compared to c = 0.02 ± 0.01 
for 0.3 < 2; < 1.1, and no significant offset in the Hubble 
diagram. 

There could however still be unresolved NICMOS is- 
sues. For example the NICMOS SN photometry depends 
on extrapolating the non-linearity correction to low flux 
levels. We have a program (HST GO-11799) to obtain a 
calibration of NICMOS at low flux levels. The photom- 
etry of the SNe observed with NICMOS will be revised 
once this program is completed. The new data presented 
in this paper also allow us another route to check this 
color discrepancy with IR data independent of the NIC- 
MOS calibration. 

7.2. Fitting Cosmology 

Following Conley et al. (2006) and K08, we adopt a 
blind analysis approach for cosmology fltting where the 
true fitted values are not revealed until the complete 
analysis framework has been settled. The blind tech- 
nique is implemented by adjusting the magnitudes of the 
SNe until they match a fiducial cosmology {Qm = 0.25, 
w = —1). This procedure leaves the residuals only 
slightly changed, so that the performance of the analy- 
sis framework can be studied. The best fitted cosmology 
with statistical errors is obtained through an iterative 
X'^-minimization of 



Astat 



= E 



[lJ,B{a,/3,M) - n{z;nM,^w,w)] 



(4) 
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Table 6 

Light curve fit results using SALT2. The SALT2 parameters are MJDmax, "Jb, ^o, and c. Here, the 
light curve shape, xi, and color, c, corrected peak B-band magnitude is obtained from equation (3) as 
m'^" = mg*™ + a ■ XI — P ■ c, where the values a = 0.12 and /3 = 2.51 have been used. For comparison 
we also present the SALT stretch, s, which has been derived from xi using the relation provided in Guy 

et al. (2007). 



SN 


MJDmax 






c 


XI 


s 


2001CW 
2001gn 
2001go 
2001gq 
2001gy 
200 Ihb 


52066.60 ±0.86 
52030.33 ±0.16 
52010.99 ±0.35 
52028.81 ±0.42 
52031.66 ±0.25 
52038.02 ± 0.97 


24.88 ± 0.42 
25.16 ±0.23 
23.13 ±0.11 
23.62 ± 0.19 
22.97 ± 0.07 
24.84 ±0.13 


24.71 ±0.13 
25.37 ±0.05 
23.08 ±0.07 
23.79 ± 0.03 
23.05 ± 0.03 
24.79 ± 0.03 


-0.10 ±0.36 
0.07 ±0.05 

-0.11 ±0.07 
0.06 ±0.04 

-0.03 ±0.03 
0.00 ± 0.04 


0.140 ±0.844 

0.648 ± 0.649 
-1.149 ±0.259 

0.459 ±0.287 
-0.312 ±0.265 

1.292 ±0.410 


0.993 ± 0.077 
1.040 ±0.061 
0.881 ± 0.022 
1.022 ± 0.027 
0.952 ± 0.024 
1.101 ±0.039 



where, 

4 = Vms+a''V,r+P^Vc+2aVmn.x^-2l3Vms,c-2apV,,,c 

is the propagated error from the covariance matrix, V , of 
the hght curve fits, witli a and [3 being the xi and color 
correction coefficients of equation (3). Uncertainties due 
to host galaxy peculiar velocities of 300 km/s and uncer- 
tainties from Galactic extinction corrections and gravi- 
tational lensing as described in 7.3 are included in (Text- 
A floating dispersion term, (Tgys, which contains poten- 
tial sample-dependent systematic errors that have not 
been accounted for and the observed intrinsic SN la dis- 
persion, is also added. The value of agys is obtained by 
setting the reduced to unity for each sample. Com- 
puting a separate agys for each sample prevents samples 
with poorer-quality data from increasing the errors of 
the whole sample. This approach does however still as- 
sume that all SNe within a sample are measured with 
roughly the same accuracy. If this is not the case there 
is a risk in degrading the constraints from the sample 
by down weighting the best measured SNe. It should 
also be pointed out that the fitted values of cTgys will be 
less certain for small samples and can therefore deviate 
significantly from the average established by the larger 
samples (in particular, the six high- 2: SNe presented in 
this work are consistent with crgys = 0), as are three other 
samples. 

A number of systematic errors are also being consid- 
ered for the full cosmology analysis. These are taken into 
account by constructing a covariance matrix for the en- 
tire sample which will be described below in 7.3. The 
terms in the denominator of equation (4) are then added 
along the diagonal of this covariance matrix. 

Following K08, we carry out an iterative minimiza- 
tion with 3a outlier rejection. Each sample is fit for its 
own absolute magnitude by minimizing the sum of the 
absolute residuals from its Hubble line (rather than the 
sum of the squared residuals). The line is then used for 
outlier rejection. This approach was investigated in de- 
tail in K08, and it was shown with simulations that the 
technique is robust and that the results are unaltered 
from the Gaussian case in the absence of contamination 
and that in the presence of a contaminating contribution, 
its impact is reduced. Table 7 summarizes the effect of 
the outlier cut on each sample. We also note that the 
residuals have a similar distribution to a Gaussian in that 
^ 5% of the sample is outside of 2a. 

Figure 7 shows the individual residuals and pulls from 



the best fit cosmology together with the fitted SALT2 
colors for the different samples. The photometric quality 
is illustrated by the last column in the figure showing 
the color uncertainty. It is notable how the photometric 
quality on the high redshift end has improved from the 
K08 analysis. This is due to the extended rest-frame 
range of the SALT2 model compared to SALT. 

Figure 8 shows the diagnostics used for studying the 
consistency between the different samples. The left panel 
shows the fitted asyg values for each sample together with 
the RMS around the best fitted cosmology. The intrinsic 
dispersion associated with all SNe can be determined as 
the median of agys as long as the majority of the sam- 
ples are not dominated by observer-dependent uncertain- 
ties that have not been accounted for. The median (jgys 
for this analysis is 0.15 mag, indicated by the leftmost 
dashed vertical line in the figure. The two mid-panels 
show the tensions for the individual samples, by com- 
paring the average residuals from the best-fit cosmology. 
The two panels show the tensions without and with sys- 
tematic errors (described in 7.3) being considered. Most 
samples fall within la and no sample exceeds 2a. The 
right panel shows the tension of the slopes of the resid- 
uals as a function of redshift. This test may not be very 
meaningful for sparsely sampled data sets, but could re- 
veal possible Malmquist bias for large data sets. 

The SNe introduced in this work show no significant 
tension in any of the panels. The Hubble residuals for 
these are also presented in Figure 9. Here, the individual 
SNe are consistent with the best fit cosmology. 

All tables and figures, including the complete covari- 
ance matrix for the sample, are available in electronic 
format on the Union webpage'^^. We also provide a Cos- 
moMC module for including this supernova compilation 
with other datasets. 

7.3. System.atic errors 

The K08 analysis split systematic errors into two cat- 
egories: the first type affects each SN sample indepen- 
dently, the second type affects SNe at similar redshifts. 
Malmquist bias and uncertainty in the colors of Vega 
are examples of the first and second type, respectively. 
Typical numbers were derived for both of these types 
of systematics, and they were included as covariances^^ 

http : //supernova . Ibl . gov/Union/ 

Note that adding a covariance is equivalent to minimizing over 
a nuisance parameter that has a Gaussian prior around zero; the 
discussion in K08 is in terms of these nuisance parameters. This is 
further discussed in Appendix C. 
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Figure 7. Individual diagrams and distributions for tiie different data sets. From left to right: a) Hubble diagrams for the various samples; 
b) binned magnitude residuals from the best fit cosmology (bin- width: Az = 0.01); c) unbinned magnitude residuals from the best fit; 
d) histogram of the residuals from the best fit; c) pull of individual SNe as a function of redshift; f) histogram of pulls; g) SN color as a 
function of redshift; h) uncertainty of the color measurement as an illustration of the photometric quality of the data. 
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Figure 8. Diagnostics plot for the individual data sets. Prom left to right: Systematic dispersion (filled circles) and RMS around the best 
fit model (open circles); The mean, sample averaged, deviation from the best fit model; The slope of the Hubble-residual (in magnitudes) 
versus redshift, d/iresidual/<i2. Note that the errors on the systematic dispersion are the statistical errorbars and do not include possible 
systematic effects such as misestimating photometry errors. 
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Table 7 

Statistics of each sample with no outlier rejection or Scr outlier rejection (used in this paper). 
Here, crsys has the same meaning as in equation (4). Both asys and the RMS are also plotted 
for each sample in Figure 8. Although each sample is independently fit for its Csys and RMS, a 
global a and l3 are always used. This explains the minor shifts in parameters for samples where 
no supernovae are cut. A 2<t cut removes 34 more supernovae, so going from 3 to 2<t is 
consist out w'lfh Gaussian residuals. 



No Outlier Cut o"cut = 3 

Set SNe (Tsys(68%) RMS (68%) SNe <T^ys{68%) RMS (68%) 



Hamuy et al. (1996) 
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Figure 9. Upper panel: Hubble diagram for the Union2 compi- 
lation. The solid line represents the best fitted cosmology for a 
flat Universe including the CMB and BAO constraints discussed 
in the text. The different colors have the same interpretation as 
in Figures 7 and 8. Lower panel: Hubble diagram residuals where 
the best flttcd cosmology has been subtracted from the light curve 
shape and color corrected peak magnitudes. The gray points show 
the residuals for individual SNe, while the black points show the 
binned values in redshifts bins of 0.05 for z < 1.0 and 0.2 for 
2 > 1.0. The orange points show the previously unpublished SNe 
introduced in this work. The dashed lines show the expected Hub- 
ble diagram residuals for cosmological models with ui ± 0.1 from 
the best fitted value. 

between SNe. Each sample received a common covari- 
ance, and all of the high-redshift {z > 0.2) SNe shared 
an additional common covariance. 

Other analyses (Astier et al. (2006), Wood-Vasey et al. 
(2007), KS09) have estimated the effect on w for each sys- 
tematic error and summed these in quadrature. However, 
Kim & Miquel (2006) show that parameterizing system- 
atic errors (such as uncertain zeropoints) with nuisance 
parameters is a more apropriate approach and gives bet- 
ter cosmological constraints. For this analysis, all con- 
tributing factors, described below, were translated to nui- 
sance parameters, which were then incorporated into a 
covariance matrix for the distances of the individual SNe. 
Appendix C contains the details of converting nuisance 
parameters to a covariance matrix. 

7.3.1. Zero-point Uncertainties 

In order to correctly propagate calibration uncertain- 
ties, we computed numerically the effect of each photo- 
metric passband on the distance modulus. For each SN, 
the photometry from each band was shifted in turn by 
0.01 magnitudes and then refit for /i. We then computed 

the change in distance modulus, giving '^^^^^p'^'* for each 
band. A list of zero-point uncertainties is given in Ta- 
ble 8. For two SNe, i and j, with calibrated photometry 
obtained in the same photometric system, the zero-point 
uncertainty, azp of that system was propagated into the 
covariance matrix element Uij as ^j^py^j^pyCzp accord- 



Table 8 

Assumed zero-point uncertainties for SNe in the Union2 
compilation. 



Source 


Band 




Reference 


HST 


WFPC2 


0.02 


Heyer et al. (2004) 




ACS 


0.03 


Bohlin (2007) 




NICMOS 


0.03 


Thattc ct al. (2009) 


SNLS 


g, r, I 


0.01 


Astier et al. (2006) 




z 


0.03 




ESSENCE 


R, I 


0.014 


Wood-Vasey et al. (2007) 


SDSS 


u 


0.014 


Kessler et al. (2009) 




9, r, i 


0.009 




z 


0.010 




'I'his paper 


R, I 


0.03 






J 


0.02 




Other 


(7-band 


0.04 


Hicken et al. (2009a) 




Other Band 


0.02 


Hicken et al. (2009a) 



ing to Appendix C.2. 

This procedure is a more efficient way of including zero- 
point uncertainties than including a common magnitude 
covariance (multiplicative in flux space) when performing 
all of the light curve fits. In testing, both of these meth- 
ods gave results that agreed at the couple of a percent 
level. Our method has the advantage that the zero-point 
errors can be adjusted without rerunning the light curve 
fits. 

Zero-point uncertainties are one of the largest system- 
atic errors (see Table 9). However, we should note that 
this number is based on a heterogeneous assessment of 
errors from different datasets (Table 8); the accuracy will 
vary. 

7.3.2. Vega 

Astier et al. (2006) estimated the broadband Vega 
magnitude system uncertainty to be within 1 % by com- 
paring spectroscopy from D. S. Hayes, L. E. Pasinetti, 
& A. G. D. Philip (1985) and Bohlin & GiUiland (2004). 
In their analysis, only the uncertainties of Vega colors 
had implications for cosmological measurements, which 
they chose to include by adopting a flux uncertainty lin- 
ear in wavelength that would offset the Vega B — R color 
by 0.01. The uncertainty of Vega is the single largest 
source of systematic error when estimating w, as shown 
in Table 9, suggesting that a better-understood reference 
would allow for a significant reduction in systematic er- 
rors. 

KS09, and recently SNLS (Regnault et al. 2009), chose 
BD+17°4708 as their primary reference star. This star 
has the advantage of having a well-known SED, measured 
Landolt magnitudes (in contrast to Vega) and colors that 
are close to the average colors of the Landolt standards 
(in contrast to Vega which is much bluer). KS09 studied 
the implications of switching between BD+17°4708 and 
Vega and found zeropoints consistent to 1%. 

Given this small difference between using BD-|-17°4708 
and Vega, we have chosen, for this work, to continue us- 
ing Vega as our primary reference star. To account for 
the uncertainty of the magnitude of Vega on the Lan- 
dolt system, we have assumed a correlated uncertainty 
of 0.01 mag for all photometry with a rest-frame wave- 
length in each of six wavelength intervals defined by the 
following wavelength boundaries: 2900A, 4000A, 5000A, 
6000A, 7000A, lOOOOA, 16000A. 
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7.3.3. Rest-frame U -Band 

SNe la are known to show increasing spectroscopic and 
photometric diversity for wavelengths shorter than the 
rest-frame i?-band. Part of this could perhaps be ex- 
plained by differences in progenitor metallicity (Hoeflich 
et al. 1998; Lentz et al. 2000), but the spectral variations 
in the rest-frame UV (Ellis et al. 2008) are larger than 
predicted by existing models. 

As discussed in Section 6, KS09 studied how well the 
SALT2 model describes the rest-frame [/-band by first 
running SALT2 with the rest frame J7-band excluded. 
Using these fits, they then generated a model for the 
rest frame J7-band and binned the magnitude residuals 
from the actual rest-frame U -band data as a function of 
phase. For the SDSS and SNLS datasets, the residuals 
around the time of maximum are ~ 3%. For SNLS, this 
is not surprising, as the SNLS data was used to train the 
SALT2 model. In this analysis, we use the SDSS sample 
as a validation set, and include a correlated 0.03 mag- 
nitude uncertainty for all photometric bands bluer than 
rest-frame 3500A. 

We note that the HST and low-rcdshift datasets are 
less useful for assessing the size of rest-frame U -band un- 
certainty. For the HST data, the light curves are poorly 
constrained without the rest-frame J7-band. In the case 
of the nearby sample, the rest-frame J7-band overlaps 
with the observed [/-band for which accurate photome- 
try is generally difficult to obtain (any potential prob- 
lems with the nearby [/-band will not impact the light 
curve fits much, as the low- .2 fits are typically very well 
constrained with the remaining bands). 

7.3.4. Malmquist bias 

K08 added a 0.02 magnitude covariance for each sam- 
ple representing Malmquist bias uncertainty. More re- 
cently, KS09 completed a thorough simulation of selec- 
tion effects for each of the samples in their analyses. 
They find a 0.024 change in w when making a correction 
for selection effects. As the 0.02 magnitude covariance 
yields a quite similar 0.026 error on w, and conducting a 
full simulation is beyond the scope of this work, we reuse 
the covariance from K08. 

7.3.5. Gravitational lensing 

The effects from gravitational lensing on the Hubble 
diagram have been discussed in detail in the literature 
(Sasaki 1987; Linder 1988; Bergstrom et al. 2000; Aman- 
uUah et al. 2003; Holz & Linder 2005). Gravitational 
lensing only affects the high redshift end of the data that 
is currently available, and potential bias on the cosmolog- 
ical parameters from the analysis carried out here due to 
the asymmetry of the lensing probability density func- 
tion is expected to be negligible. We adopt the K08 
approach of only treating gravitational lensing as a sta- 
tistical uncertainty by adding a value of 0.093z (Holz & 
Linder 2005) in quadrature to (Text in equation (4). This 
is a conservative approach with respect to the values pre- 
sented by Jonsson et al. (2006), where they attempt to 
measure the lensing of individual SNe by determining the 
mass distribution along the line of sight. A very impor- 
tant conclusion of their work is that there is no evidence 
for selection effects due to lensing of the high-redshift 
SNe. 



7.3.6. Light curve model 

We have studied any potential bias that could arise 
from poor light curve sampling, by carrying out a similar 
analysis to K08, updated for SALT2. We use nine BVR 
AQUAA templates (Strovink 2007) constructed from ob- 
servations of very well-observed nearby supernovae. Each 
set of BVR templates is combined with a SALT2 f/-band 
template generated for that supernova, as there were 
insufficient observations in the J7-band to construct an 
AQUAA template. 

Mock data sets are then sampled from these templates 
with the same rest frame dates and signal-to-noise ra- 
tios of the real SNe in our sample. The mock sets are 
then fitted with SALT2 and the offset between AQUAA 
corrected magnitude and the corresponding SALT2 fit- 
ted value is investigated as a function of the fitted phase 
of each supernova's first data point (the phase is with 
respect to U-band maximum) . K08 looked at other pos- 
sible biases, but first phase was the only significant one 
found. The test is carried out for each of the nine tem- 
plate SNe. 

For SNe with a first phase at i?-band maximum, the 
average bias is close to zero, with an RMS of about 0.03. 
For SNe with a first phase at six days past B-band maxi- 
mum, the average bias is still close to zero, but the RMS 
has increased to about 0.08. Of course, our nine SN tem- 
plates might not be a representative sample, but these re- 
sults are encouraging, since they both suggest that there 
is no significant bias and indirectly validate the SALT2 
performance with respect to the AQUAA templates. 

We use a first phase cut of six days, but we conserva- 
tively give each SN that has first phase greater than zero 
a 0.03 magnitude covariance. Note that SALT2 does not 
stretch its definition of phase with light curve width. 

7.3.7. Contamination 

As already mentioned, we perform an iterative min- 
imization with a 3a outlier rejection before fitting cos- 
mology. In K08 we showed that this technique greatly re- 
duces the impact of potential contamination, while main- 
taining roughly Gaussian statistics. We carried out a 
Monte-Garlo study showing that the effect of contami- 
nation on any individual sample is limited to less than 
0.015 magnitudes. This is under the assumption that 
the dispersion of the contaminating distribution is of the 
same order, or greater than, the dispersion of SNe la and 
that the contamination is less than 30 %. We include a 
0.015 magnitude uncertainty, correlated for each sample, 
to account for possible contamination. 

7.3.8. Minimum redshift 

In order to test for possible effects from using a given 
minimum redshift cut, we started by constructing a new 
sample with no minimum redshift. Using this sample, we 
performed fits which allowed the absolute magnitude to 
vary independently below and above a dividing redshift 
in the range 0.01 to 0.03. This procedure should test for 
a Hubble bubble or significantly correlated peculiar ve- 
locities. The extra degree of freedom allowed by this step 
in Mb improved the by < 1 regardless of the divid- 
ing redshift and the inclusion of systematic errors. This 
confirms the results of Gonley et al. (2007) for SALT2. 
We conclude that there is no statistically significant dif- 
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ference between minimum redshifts and use the value of 
0.015, as was used in the K08 analysis. 

7.3.9. Galactic extinction 

All light curve photometry is corrected for Galactic 
extinction using the extinction law from Cardelli et al. 
(1989), assuming Ry = 3.1, together with the dust maps 
from Schlegel et al. (1998). 

In the same procedure as with calibration uncertain- 
ties, we increased the Galactic E{B — V) by 0.01 for 

each supernova and repeated the fit, giving j^i^^b-v)) • 
A 10% statistical and 16% systematic error was assumed 
for the Galactic extinction of each supernova (Schlegel 
et al. 1998). 

7.3.10. Intergalactic extinction 

Dimming of SNe la by hypothetical intergalactic gray 
dust has been suggested by Aguirre (1999) as an alter- 
native to dark energy to explain the SN results (Goobar 
et al. 2002). This potential dimming was however con- 
strained by studying the colors of high-z quasars (Mort- 
sell & Goobar 2003; Ostman & Mortsell 2005) and by ob- 
servations SNe la in the rest frame I-band (Nobili et al. 
2005, 2009). 

Another possible extinction systematic comes from the 
dust in galaxy halos that are along the line of sight. 
Menard et al. (2009b) used distant quasars to detect and 
measure extinction in galactic halos at z ^ 0.3. They 
find an average Ry for their galaxies of 3.9 ± 2.6. Us- 
ing their observed Ayir), we find an average rest frame 
V-band extinction of 0.004 magnitudes per intersected 
halo, assuming Ry =3.1. At redshift 0.5, an average of 
three halos have been intercepted. At redshift 1.0, the 
average is seven. 

There are three mitigating factors. One is that expan- 
sion redshifts photons between the supernova and the 
intervening galaxy. The COM law decreases with wave- 
length (in the relevant wavelength range), so less light 
is absorbed. Menard ct al. (2008) finds that pdust k 
(1 -|- 2;)^^'^, which we use to scale the extinction. Finally, 
most of the extinction is corrected by color correction. 
The exact amount corrected depends on the redshift and 
the filters used in the observations, but is around two 
thirds. 

We find an error on w of 0.008 due to this extinc- 
tion, significantly lower than the value of 0.024 derived 
by Menard et al. (2009a). However, they used Ry = 3.9, 
rather than 3.1; the fraction of extinction that is cor- 
rected by the color correction will decrease with Ry. We 
also numerically sum the CCM laws, rather than using 
an analytic approximation. Since we know the exact red- 
shift and filters used in each observation, we can exactly 
calculate the amount of extinction already handled by 
the color correction (using our ^ values), without ap- 
proximation. 

7.3.11. Shape and Color Correction 

The most uncertain contribution to the dimming of 
SNe la is host galaxy extinction. Several studies of SN la 
colors (Guy et al. 2005, 2007; Wang et al. 2008; Nobih 
& Goobar 2008, and references therein) indicate that the 
observed SN la reddening does not match the Galactic 



CCM extinction law with Ry = 3.1. A stronger wave- 
length dependence has been found in the optical in most 
cases, and it remains unclear if CCM models with any 
value of Ry can be used to describe the data accurately. 
It is possible that the observed steep reddening originates 
from a mixture of local effects and host galaxy extinc- 
tion. Local effects could be intrinsic SN variations, but 
also multiple scattering on dust in the circumstellar envi- 
ronment has been suggested (Wang 2005; Goobar 2008). 
This model is potentially supported by detection (Patat 
et al. 2007; Simon et al. 2009; Blondin et al. 2009) of 
circumstellar material but also by color excess measure- 
ments for two of the best observed reddened SN la (Fo- 
latelli et al. 2010) being consistent with the expected 
extinction from circumstellar material. 

The SALT2 method approaches the lack of a consistent 
understanding of SN la reddening by adopting a purely 
empirical approach. For SALT2, the SN la luminosity is 
standardized by assuming that the standardization is lin- 
ear in both xi and c as described in equation (3), where /3 
is the empirically determined correction coeffecient that 
accounts for all linear relations between color and ob- 
served peak magnitude. For example if the only source 
for such SN la reddening originated from CCM extinc- 
tion then (3 is identically equal to Ry + 1. We test this 
approach and propagate relevant systematic unccirtain- 
ties by dividing the full sample into smaller sets and car- 
rying out independent fits for the Xi and c correction 
coefficients, a and /9, as shown in Table 11. 

When subdividing into redshift bins, we find that the 
values of a and /3 for the full sample are consistent with 
values for the three first redshift bins. It is encouraging 
to see consistency between the global values fit for the 
full dataset and the values in the best-understood red- 
shift range. Beyond this important test, we also note 
that the value of (3 in the redshift range 0.5 to 1 is signif- 
icantly lower than the other values, while the value for 
2: > 1 is higher than the global value, but is poorly mea- 
sured. The trend is similar to what was seen in KS09, 
but we use different binning. This behavior is inconsis- 
tent with a monotonic drift in redshift, so we consider 
other explanations for these results. That conclusion is 
also supported by the observation that samples at similar 
redshifts (e.g. Miknaitis et al. (2007) and Astier et al. 
(2006)) can have very different values of (3 when fitted 
independently. The value of a is consistent across red- 
shift ranges, except at z > 1, where many light curves are 
so poorly sampled that it may not be possible to assign 
reasonable xi errors. 

When subdividing into the four largest data sources 
(the lower half of Table 11), we find values of a and /3 
generally consistent with the global values, with the ex- 
ception of a lower value of for the SNLS SNe (Astier 
ct al. 2006). In general, ignoring or underestimating the 
errors in c or xi will decrease the associated correction 
coefficients, /? and a, as investigated in K08 and KS09 
and this may be relevant here. Specifically, two potential 
sources of problems are an incomplete understanding of 
calibration and underestimated SN model variations, ei- 
ther of which could affect these fits. If the SNLS SNe are 
physically different, and they are allowed their own /?, 
then w shifts by 0.02. Alternatively, as one f3 is used for 
the global sample, the possibility that that (3 is biased 
from the true global value must be considered. Selecting 
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the global value of j3 from any of the other large sam- 
ples shifts w by less than 0.02. We have accounted for 
this systematic by assigning each sample a 0.02 magni- 
tude covariancc (giving a 0.03 error on w), which avoids 
the problem of handling an error on elements of the co- 
variance matrix. To study these details further, we look 
forward to more data for z > 0.5, with improved calibra- 
tion and light curve models. 

We also perform one additional sanity check by sub- 
dividing the data by xi and c. There is evidence for 
two populations of normal SN la, divided by light curve 
width (sec K08, and references therein). Star- forming 
galaxies tend to host the population with broader light 
curves, while SN hosted by passive galaxies tend to have 
narrower light curves. As described below, we derive 
consistent cosmology for these subdivisions as well. 

We subdivide^^ the full sample into two roughly equal 
subsamples, split first by color and then by xi. In total, 
this makes four subsamples. We find that the cosmology 
is close for all subsamples (as can be seen in Table 11) so 
the difference from these subdivisions does not contribute 
significantly to the systematic error on w. 

It is interesting to note that a is substantially different 
for the two samples split by light curve width. Likewise, 
13 is substantially different for the two samples split by 
color. This might suggest that the relationships between 
color and brightness and light curve width and brightness 
are more complex than a simple linear relationship, or it 
could be that the errors on xi and c are not perfectly 
understood. We also find that /? is higher for the redder 
SNc la which is similar to the results from Conley et al. 
(2008) based on comparisons oi U — B colors to i? — 
V, after correcting for the effect of stretch on the U- 
band. At the same time it should be pointed out that 
evidence of low Ry values have also been found for a few 
well-studied, and significantly reddened, nearby SNe la 
Folatelh et al. (2010). 

7.3.12. Summary of systematic errors 

The effect of these systematic errors on w is given in 
Table 9. The improvement in cosmology constraints over 
the simple quadrature sum is also shown. Zeropoint and 
Vega calibration dominate the systematics budget, but 
understanding the color variations of SNe is also impor- 
tant. The benefit from making a Malmquist bias correc- 
tion can be seen; by doing so, KS09 reduce this system- 
atic error by a factor of two. 

8. RESULTS AND DISCUSSION 

In the cosmology analysis presented here, the statisti- 
cal errors on ^Im have decreased by a significant 24% over 

Subdividing by xi or c must be done carefully. When there 
are errors in both the dependent and independent variables (in this 
case, magnitude and xi or c), the true values of the independent 
variables must be explicitly solved for as part of the fit. Other- 
wise, the subdividing will be biased. For example, suppose that 
a supernova has a color tfiat is poorly measured, and an uncor- 
rected magnitude that is well-measured. If this supernova is faint 
and blue, then a fit for the true color will give a redder color. A 
color cut will place this supernova in the blue category, when the 
supernova is actually more likely to be red. As mentioned in K08, 
whenever one fits for a and /3, the true values of a;i and c are only 
implicitly solved for; equation (5) is derived by analytically mini- 
mizing over the true xi and c. K08 provides the equation with the 
true values made explicit, we also include a discussion in Appendix 
C. 



Table 9 

Effect on w errorbar (including BAO and CMB 

constraints) for each of the systematic errors 
included. The proper way to sum systematic 
errors is to include each error in a covariance 
matrix. 



Source 



Error on w 



Zero point 0.037 

Vega 0.042 

Galactic Extinction Normalization 0.012 

Rest-Frame [/-Band 0.010 

Contamination 0.021 

Malmquist Bias 0.026 

Intergalactic Extinction 0.012 

Light curve Shape 0.009 

Color Correction 0.026 



Quadrature Sum (not used) 
SuiiiiiiocI ill CoA'aiiaiiee Matiix 



0.073 
0.063 



the K08 Union analysis, while the estimated systematic 
errors have only improved by 13%. When combining the 
SN results with BAO and CMB constraints, statistical 
errors on w have improved by 16% over K08, though the 
quoted systematic errors have increased 7%. Figure 10 
shows a comparison between the constraints from K08 
and this compilation in the (Om — w) plane. Even with 
some improvement on the understanding of systematic 
errors, it is clear that the dataset is dominated by sys- 
tematic error (at least at low to mid-z). 

The best fit cosmological parameters for the compi- 
lation are presented in Table 10 with constraints from 
CMB and BAO. The confidence regions in the (JIm, ^^a) 
and (r^M , w) planes for the last fit in the table are shown 
in Figures 11 and 12 respectively. 

For the CMB data we implement the constraints from 
the 7 year data release of the Wilkinson Microwave 
Anisotropy Probe (WMAP) as outlined in Komatsu et al. 
(2010). We take their results on z^, (the redshift of last 
scattering), Z^(z^), and R{Zi,), updating the central val- 
ues for the cosmological model being considered. Here, 
lA{Zi,) is given by 



Ia{z^) = (1 + -2*) 



rsizi,) 

where Da is the angular distance to z*, while 



{1 + z^)Da{z^) . 



Percival et al. (2010) measures the position of the BAO 
peak from the SDSS DR7 and 2dFGRS data, constrain- 
ing = rs{zd)/Dvi0.275) to 0.1390 ± 0.0037, where 
fz{zd) is the comoving sound horizon and Dv{z) = 

[{l + zrD\cz/H{z)f'. 

For the SNe + BAO fit in Table 10, we add an Hq mea- 
surement of 74.2±3.6 km/s/Mpc from Riess et al. (2009), 
creating a constraint without the CMB that is therefore 
largely independent of the high-redshift behavior of dark 
energy (as long as the dark energy density contribution 
is negligible in the early universe). Note that the Hq 
constraint relies on most of the nearby supernovae used 
in this compilation. However, the effect on w through 
Hq from these supernovae is several times smaller than 
the effect through the Hubble diagram. Alternatively, 




Table 10 

Fit results on cosmological parameters Qm, ^ and Ofe. The parameter values are followed by their statistical 
(first column) and statistical and systematic (second column) uncertainties. 



Fit 




flAf w/ Sys 




flk w/ Sys 


w 


w w/ Sys 



SNe 0.2701°;q21 0.274l°-°3? (fixed) (fixed) -1 (fixed) -1 (fixed) 

SNe+BAO+i^o 0.3091°;°^^ O.Sietoiogg (fixed) (fixed) -l-114to.??2 -'^■'^^'^-oaIo 

SNe+CMB 0.268tJ|!;E;i? 0.269tU;!522 (S^ed) (fixed) -0.997t|5;!5^^ -0.999tim 

SNe+BAO+CMB 0.277to:o" 0.279t°°ll (fixed) (fixed) -'^■OO^to°o54 -^■^^'''-omII 

SNe+BAO+CMB 0.278too" 0.28lt°-g}^ -0-004l°°°^ -0.004tooo? -1 (fixed) -1 (fixed) 

SNe+BAO+CMB 0.28lto;EJ}^ 0.28lto;«i« -0-005tiZ -O-OOetlZ' "l-OSe^-.^fg -1.0351^9? 



adding a CMB constraint on O^/i^ of 0.1338 ± 0.0058 
from the WMAP7 webpage"^^ allows us to create a con- 
straint that is independent of Hq. This final result for 
SNe +BAO-I-CMB does not improve significantly if the 
current Hq constraint is added. 

8.1. Time variation of the dark energy equation of state 

The constraints shown in Figure 12 were obtained as- 
suming that the dark energy equation of state (EOS) is 
redshift independent. SNe la are useful for constraining 
a redshift dependent w{z) since, unlike e.g. CMB, their 
measured distances at a given redshift are independent 
of the behavior of dark energy at higher redshifts. A 
common method to parameterize w{z) is 

W{Z) ^ Wo+ Wa— 

1 + z 

where a cosmological constant is described by {wo,Wa) = 
(— 1, 0). It can be shown (Linder 2003) that this parame- 
terization provides an excellent approximation to a wide 
variety of dark energy models. The constraints from the 
current SN data together with the CMB and BAO data 
are presented in Figure 13. The flattening of the contours 

^® http : //lambda . gsf c . nasa . gov/product/map/current/params/wc 



in this diagram at wq + Wa — comes from the implicit 
constraint of matter domination in the early Universe 
imposed by the CMB and BAO data. Only modest con- 
straints can currently be placed on Wa- 
it can be illuminating to study w{z) in redshift bins, 
where w is assumed constant in each bin. This method 
has the advantage that w{z) can be studied without as- 
suming a specific form for the relation (Huterer & Stark- 
man 2003). We carry out the analysis following Kowalski 
et al. (2008) and fit a constant w in each bin, while the re- 
maining cosmological parameters are fit globally for the 
entire redshift range. Figure 14 shows three such mod- 
els for the combined constraints from SNe, BAO, CMB, 
and Hq measurements where we assume a flat Universe. 
In these scenarios, the Hq measurement does not con- 
tribute much, but due to its small improvement on the 
CMB constraints it gives a small (^ 10%) improvement 
on the errorbar of the highest redshift bin. 

The left panel shows constraints on w for three bins. 
The first bin (0 < z < 0.5) shows a well-constrained w. 
The middle bin (0.5 < z < 1) shows a poorly-constrained 
w, though one that is distinct from — oo (which would 
drive p to above z = 0.5, resulting in a matter-only 
universe) at high confidence, indicating the detection of 
, some kind of dark energy in this redshift range. For 

(im_sz_lens_winap7 . cim ° 
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Figure 11. 68.3%, 95.4%, and 99.7% confidence regions in tlie {Qm,^Ia) plane from SNe combined witli the constraints from BAO and 
CMB both without (left panel) and with (right panel) systematic errors. Cosmological constant dark energy (w = —1) has been assumed. 



z > 1, there is little constraint on w, and only a weak 
constraint on the existence of dark energy. 

The middle panel shows the effect of dividing the high- 
est redshift bin. The constraints on w for z > 1 get much 
weaker, showing that most of the (weak) constraint on 
the highest bin in the left panel comes from a combina- 
tion of the CMB with the well-constrained low-redshift 
supernova data. Current supernovae at z > 1 offer no 
real constraint on w(z > 1). Providing a significant con- 
straint at these redshifts requires significantly better su- 
pernova measurements. As in the left panel, w in the 
highest redshift bin is constrained to be less than zero 
by the requirement from BAO and CMB constraints that 
the early universe have a matter-dominated epoch. 

The right panel shows the effect of dividing the low- 
est redshift bin. While no significant change in w with 
redshift is detected, there is still considerable room for 
evolution in w, even at low redshift. 

Figure 15 shows dark energy density constraints, as- 
suming the same redshift binning as in Figure 14. Note 
that this is not equivalent to the left and center panels of 
Figure 14; only in the limit of an infinite number of bins 
do binned p and binned w give the same model. Dark 
energy can be detected at high significance in the middle 
bin (redshift 0.5 to 1), but there is only weak evidence 
for dark energy above redshift 1 (left panel). When the 
bin above redshift 1 is split at a redshift greater than the 
supernova sample (right panel), it can be seen that the 
current small sample of supernovae cannot constrain the 
existence of dark energy above redshift 1. 



8.2. SNe with ground-based near-IR data 

Obtaining near-IR data of z > 1 SNe la, whether from 
space or from the ground, is critical for constraining the 
SALT2 color parameter, c. Without the near-IR data, 
the uncertainty in this parameter for 2001hb and 2001gn, 
both beyond z = 1, increases by a factor of two. Precise 
measurements of c are important, since uncertainties in 
c are infiated by (3 w 2.5 and tend to dominate the error 
budget when the corrected peak i?-brightness of SNe la 
are calculated. 

Both 2001hb and 2001gn were observed with ground- 
based near-IR instruments. The operational challenges 
associated in obtaining these data are significant. Long 
exposure times (ten hours or more taken within a few 
days) in excellent observing conditions are necessary. 
Even with queue mode scheduling, these observations are 
just feasible. Despite the challenges, the uncertainty in 
the SALT2 color of these two SNe la is comparable to the 
uncertainty in the color of the best space-based measured 
SNe la at z > 1. 

The ground based near-IR data also allow us to search 
for systematic offsets with near-IR data taken from space. 
For z > 1.1 SNe la observed with NICMOS, the average 
SALT2 c value is c = 0.06 ± 0.03 mag. By compari- 
son, the weighted average color of the three SNe la at 
z 1.1 with ground-based near-IR data (2001hb and 
2001gn from this work, together with 1999fk from Tonry 
et al. (2003)) that pass the light curve cuts is, 0.01 ±0.07. 
Neither the ground-based or space-based measurements 
show any Hubble diagram offset, (A/it = 0.03 ± 0.10 and 
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Figure 12. 68.3%, 95.4%, and 99.7% confidence regions of the (njv/,io) plane from SNe combined with the constraints from BAO and 
CMB both without (left panel) and with (right panel) systematic errors. Zero curvature and constant w have been assumed. 

for these SNe, we noticed that the H09 KeplerCam r- 
band photometry is generahy 0.05 ± 0.02 magnitudes 
fainter than the KS09 SDSS r-band photometry based 
on SALT2 fits including both data sets. The quoted 
uncertainty is completely dominated by the zero-point 
uncertainties, as both sets of data have high S/N. The 
offset is consistent for each supernova. This offset is the 
correct size and direction to explain the tension in Hub- 
ble residuals seen between these datasets in second panel 
of Figure 8. 

SN 2005hc and SN 2005ir were also observed by CSP 
using the Swope telescope. In these cases, the Swope 
photometry agrees with the SDSS r-band photometry. 

9. SUMMARY AND CONCLUSIONS 

We have presented the light curves of SNe la in the 
redshift range 0.511 < z < 1.12. The SNe were dis- 
covered as part of a search conducted by the Supernova 
Cosmology Project using the Subaru, CFHT and CTIO 
Blanco telescopes in 2001. The fitted light curve shapes 
and colors at maximum light are all consistent with the 
corresponding distributions from previous SN surveys. 

Following K08, we add these six SNe and other SN la 
data sets to the Union compilation. We have also im- 
proved the Union analysis in a number of respects, cre- 
ating the new Union2 compilation. The most impor- 
tant improvements are: (1) Systematic errors are directly 
computed using the effect they have on the distance mod- 
ulus (2) All SN light curves are fitted with the SALT2 
light curve fitter. 

We determine the best fit cosmology for the Union2 
compilation, and the concordance ACDM model remains 
an excellent fit. The new analysis results in a significant 
improvement over K08 in constraining w over the redshift 
interval < z < 1. Above z > 1, evidence for dark 
energy is weak. This will remain the case until there is 
much more high redshift data, with better signal-to-noise 
and wavelength coverage. 

Based, in part, on observations obtained at the ESO La 




Figure 13. 68.3%, 95.4%, and 99.7% confidence regions of the 
(wo, Wa) plane from SNe combined with the constraints from BAO 
and CMB both with (solid contours) and without (shaded con- 
tours) systematic errors. Zero curvature has been assumed. Points 
above the dotted line (uiq -I- Wa = 0) violate early matter domina- 
tion and are implicitly disfavored in this analysis by the CMB and 
BAO data. 

—0.01 ± 0.06, respectively), from the best fit cosmology. 
These results include a 0.1 magnitude dispersion in color, 
and the fitted systematic dispersions in magnitude. Ad- 
ditional SNe la at z ^ 1.1 with ground based near-IR 
photometry will be published in Suzuki et al. (in prepa- 
ration) and we will revisit this issue in that paper, where 
we will also incorporate the results from the NICMOS 
calibration program referenced in section 7.1. 

8.3. Comparison of KeplerCam and SDSS Photometry 

H09 and KS09 share three normal supernovae in com- 
mon: 2005hc, 2005hj, and 2005ir. In comparing the data 
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Figure 14. Constraints on w{z), where w(z) is assumed to be constant in eacli rcdshift bin, are plotted at tlie 68% probability level. The 
results were obtained assuming a flat Universe for the joint data set of supernovae, BAO, CMB, and Hq, with (dark/orange) and without 
(light/yellow) SN systematics. The left panel shows three redshift bins, with the highest redshift bin keeping w constant for all z > 1. Dark 
energy is seen to exist at 2 > 1 (at least at 68 % cl) since ui does not reach an infinitely negative value, indicating its density does not go to 
zero. The middle panel splits this last bin into two, showing that the seemingly tight constraints on dark energy at 2 > 1 with current data 
depend on the combination of CMB with low-redshift data. No current probe alone can constrain the existence of dark energy at 2 > 1. 
The right panel shows the effects of w binning at low redshift. The best flt values of w go from less than —1 at 2 = 0.14 to greater than 
— 1 at 2 = 0.04. While such a steep, late time transition in w (corresponding to dw/d\na ^ 7) is unusual in physical models, it can easily 
appear due to offsets between heterogeneous data sets. We emphasize that the results are still consistent with the cosmological constant 
(dot-dashed line) at the 68% confldence level. 



SNe 

+BAO 

+CMB 

+//(, 




Figure 15. Here the dark energy density in units of the critical density today, p/pco, is assumed constant in each bin. The same binning 
as the left/center panels in Figure 14 is chosen. As can be seen in the left panel, dark energy is detected between redshift 0.5 and 1 at high 
significance, but only hints of dark energy are seen above redshift 1. When the CMB and SNe are separated, neither one alone can provide 
any evidence for dark energy at z/,1. 
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Table 11 

Subdivisions of the Union2 compilation. Values of absolute S-band magnitude, Mg (assuming 
-ffo = 70 km/s/Mpc), as well as stretch and color correction coefficients, a and /3, for several redshift ranges. 
Qm and w are shown for properly conducted xi and c cuts, chosen to give similar uncertainty of the fitted w. 
The outlier rejection is redone for each bin, so the totals may not add up to the whole sample. The constraints 

are computed including BAO and CMB data. 



Cut 


Number 


Mb 


a 





Qm 


w 



0.015 <z< 0.10 
0.100 <z< 0.25 
0.250 <z< 0.50 
0.500 <z< 1.00 
z > 1.000 



166 
74 

154 

13.3 
16 



l».OZO_Q gj^g 

i.y.d^D_|., (,2() 

-19 309+° °^^ 
-19 4:50+° °^^ 



0.112 



-19 

-19.305t°«if, 

-i9.358t;;:«i 

-19 335+0 031 
-19 315+" "^^ 
-19 299+° °^^ 
-19.325«:°?,?, 

19.31oiroit 



+0.011 
-0.011 
U.±O1_0 0jj^ 

g -, -,(.+0.013 
U.liU_0 0J3 

129+0-018 

U.1Z»_0 0jg 

-0 124+0-085 

"•-^^^-0.118 

112+'^-"^^ 
"^•-^ -^^-0.011 

n -199+0.010 

U.±ZZ_„ 0J0 

U.UZD_Q 021 
147+0.020 

'^•i^'-a.oi9 

"•^^'-0.013 
(, +0.013 
U.110_Q 013 

(, -, ,9+0-037 

140+0-017 

U.1^U_0 0iy 



9 77+0.09 
^- ' ' -0.09 

2 4q+o.l."-> 

^•"^-0.11 

1 45+0-19 

3 84+1-20 
'^•'^^-0.85 

2 96+0-10 

^•^"-0.10 

1 i., + 0.30 

2 60+0-10 
2 42+0-10 

9 70+0.10 
^•'•^-0.10 

2 50+0-1^ 
-I 79+0.18 
l-'^-0.17 



0.270 (fixed) 
0.270 (fixed) 
0.270 (fixed) 
0.270 (fixed) 
0.270 (fixed) 



o-283t;!:i- 

0.284l»;0- 

0.2781":;:- 

276+0-010 

"•^ '"-0.015 



-1.000 (fixed) 

-1.000 (fixed) 
-1.000 (fixed) 
-1.000 (fixed) 
-1.000 (fixed) 



-0 969+"-"'" 

U.,JD,J_„ 074 
3+0.058 
-0.062 

-0.996+"-«" 



c > 0.05 
c < 0.05 
XI > -0.25 
xi < -0.25 



245 

308 
302 
254 



-0.959^ 



-1.016 



-0.065 
+0.069 
0.073 



Holtzman et al. (2009) 129 

Hicken et al. (2009) 102 

Miknaitis et al. (2007) 74 

Astier et al. (2006) 71 



191+U.UU7 

U.J.^J._Q 007 



9 t;i+U:U7 



0.270 (fixed) 
0.270 (fixed) 

0.270 (fixed) 
0.270 (fixed) 



„ 9~-+U.U14 
<J.^' '-0.014 



-1.000 (fixed) 
-1.000 (fixed) 

-1.000 (fixed) 
-1.000 (fixed) 



,+U.U5U 



z > 0.015 



557 - 



-1.009 
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APPENDIX 

SPECTRA AND NOTES ON INDIVIDUAL CANDIDATES 

The spectra of the six SNe in this paper. Details of the spectra for 2001go and 2001gy can be found in Lidman et al. 
(2005), while the spectrum for 2001cw is described in Morokuma et al. (2010). 

The spectrum of each candidate is plotted twice. In the upper panel, the spectrum of the candidate is plotted in 
the observer frame and is uncorrected for host galaxy light. Telluric absorption features are marked with the symbol 
®. In the lower panel, contamination from the host (if any) is removed, and the resulting spectrum is rescaled and 
re-binned. This spectrum is plotted in black as a histogram and it is plotted in both the rest frame (lower axis) and 
the observer frame (upper axis). For comparison, low redshift SNe are plotted as a blue continuous line. 
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Figure 16. Upper panel: The error weighted and rebinned ESI spectrum of 2001gn, a probable SN la at z = 1.124. The redshift is 
derived from the [Oil] AA3727, 3729 doublet, which is clearly resolved in the ESI spectrum. Lower Panel: The rebinned host galaxy 
subtracted spectrum and SN 1990N at —7 days. 




Figure 17. Upper panel: The error weighted and rebinned ESI spectrum of 2001hb, a probable SN la at z = 1.03. Lower Panel: A 
rebinned version of the spectrum in the upper panel and SN 1989B at —7 days. The redshift of 2001hb is derived by fitting the SN to local 
SN templates as no features from the host galaxy are detected. 
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Figure 18. Upper panel: The unbinned LRIS spectrum of 2001gq, a probable SN la at z = 0.672. Lower Panel: The rebinned host 
galaxy subtracted spectrum and SN 1999bp at —2 days. 
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Table 12 

Follow-up photometry of the six SNe presented in this paper. The columns list the SN name, the 

weighted average MJD of the observation, exposure time, filter pass band, measured flux, 
instrumental ZP, and image quality in arc seconds. The ZPs are in the Landolt systems, and the 
instrumental magnitude, m, is obtained as m = — 2.51ogi() / + ZP, where / can be obtained from the 

Flux column. 



SN 


MJD 


Instrument 


Exp. 


Band 




Flux 




ZP 


I.Q. 


2001CW 


52024.43 


SuprimeCam 


3600 


i 


1.6E 


+ 00 {5.0E 


-01) 


27.30 (0.03) 


1.05 


2001CW 


52049.50 


SuprimeCam 


3240 


i 


1.3E 


+ 01 (4.6E 


-01) 


27.30 (0.03) 


1.47 


2001CW 


52079.88 


WFPC2 


2300 


F814W 


8.3E 


- 01 (4.3E 


- 02) 


23.72 (0.02) 


0.10 


200 lew 


52079.94 


WFPC2 


2400 


F850LP 


2. IE 


- 01 (4.5E 


- 02) 


22.03 (0.02) 


0.10 


2001CW 


52088.77 


WFPC2 


2300 


F814W 


5.6E 


- 01 (4.2E 


-02) 


23.72 (0.02) 


0.10 


2001CW 


52098.81 


WFPC2 


2300 


F814W 


3.4E 


-01 (4. IE 


-02) 


23.72 (0.02) 


0.10 


2001CW 


52108.91 


WFPC2 


2300 


F814W 


1.8E 


- 01 (3.7E 


-02) 


23.72 (0.02) 


0.10 


2001CW 


52118.81 


WFPC2 


4400 


F814W 


1.4E 


- 01 (2.9E 


-02) 


23.72 (0.02) 


0.10 


2001gn 


51998.47 


CFH12k 


9240 


I 


-2.3E 


+ 02 (4.0E 


+ 02) 


33.08 (0.03) 


0.80 


2001gn 


52015.50 


CFH12k 


9900 


I 


2.5E + 03 (4.5E + 02) 


33.08 (0.03) 


0.70 



Note. — This table is published in its entirety in the electronic edition of the Astrophysical 
Journal. A portion is shown here for guidance regarding its form and content. 

LIGHT CURVE DATA 

Tabulated below are light curve data for the six SNe presented in this paper. The data are listed in chronological 
order for one SN at a time, and correspond to data points plotted in Figure 6. 

All photometry is given as the measured flux, /, in the instrumental system, i.e. no color correction (S'-correction) 
to a standard system has been applied. The instrumental Vega magnitude may be calculated using the standard 
formula: m = —2.5 log / + ZP, where ZP is the associated zeropoint. For the HST data these are the updated values 
from Dolphin (2000) given on Andrew Dolphin's webpage"^''. The zeropoints have then been converted to e/s by using 
a gain of 7.12. The zeropoints for the ground-based data on the other hand are always given for the photometric 
reference image, and therefore include exposure time, gain, and for the PSF fitted data, the aperture correction. Note 
that there are correlated errors between the data points that come from the same light curve build, since the same 
host galaxy model and zeropoint has been used. This is also reflected by small offsets between different data sets in 
Figure 6. The covariances are taken into account for the light curve fitting and the full covariance matrices used will 
be available from the SCP Web site. 

The instruments used are as given below. MOSAIC II is the multichip imager on the CTIO 4 m telescope. CFH12k 
is the multichip imager on the 3.6 m CFHT telescope on Mauna Kea in Hawaii. SuprimeCam is the wide-field imager 
on the 8.2 m Subaru telescope. SuSI2 is the imager on the NTT 3.6 m telescope at ESO. FORSl is the imager on the 
8 m UT2 (Kueycn) at VLT, ESO. ISAAC and NIRI arc the near-infrared imagers on the 8 m UTl (Antu) at VLT, 
ESO and Gemini North telescopes respec:tively. Finally, WFPC2 indicates data obtained with the Planetary Camera 
CCD and and ACS is the Advanced Camera for Surveys, both on HST. 

NUISANCE PARAMETERS AS COVARIANCES 

Suppose that observations y are modeled by F{d), where some of the parameters enter into the model linearly (0^) 
while some enter non- linearly {6^). The can be written as 

X\e) = {y- (F^(0^) + He'')fV-\y - (F^(0^) + He"-)) (CI) 

where V is the covariance matrix of the observations, and H is the Jacobian matrix of the model with respect to the 
6^. Taking the derivative of the x^ and setting it to zero gives the analytic formula for the best-fit 0^ 



0^ = {H^V-'H)-'H^V-\y - FMie'')) = D{y - Fw(0^)) . 



(C2) 



The likelihoods for the 6 can be found from this restricted parameter space where = 6 . Thus, the restricted x^ 
is given by: 



X' = ((/ - HD)iy - FN{d'')) fV-\l - HD){y - F^(0^)) 



{y - FN{d'')fU-\y - F^id'')) 



with U ^ given by 



U-^ = {I - IID?'V-Hl - II D) 



y-^- - v-^ii[ii^'v-^-ir)-^-ii' v-^ 



(C3) 
(C4) 

(C5) 
(C6) 



^"^ http : //purcell . as . arizona.edu/wfpc2_calib/ 
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Table 13 

All SNe included in this compilation. The columns represent the SN name, redshift, peak 
B-band magnitude, SALT2 parameters xi and c, distance modulus, sample and failed cuts (if 

any). 



SN 


•^CMB 


^max 
B 


Xl 


c 




sample 


cut 


1993ah 


0.0285 


16.86(0.19) 


-2.26(0.93) 


u.zo 




34.61(0.23) 


1 




1993ag 


0.0500 


17.79(0.05) 


-1.09(0.24) 


0.12 


Jo. 02) 


35.95(0.17) 


1 




1993o 


0.0529 


17.60(0.05) 


-1.03(0.14) 


-0.01 


'0.02) 


36.09(0.16) 


1 




1993b 


0.0701 


18.43(0.04) 


-0.53(0.21) 


n no 


09^ 


36.71(0.16) 


1 




1992bs 


0.0627 


18.25(0.05) 


-0.27(0.23) 


0.02 


'0.02) 


36.75(0.16) 


1 




1992br 


0.0876 


19.19(0.11) 


-2.97(0.38) 


-0.04 


'0.07) 


37.50(0.19) 


1 




1992bp 


0.0786 


18.27(0.04) 


-1.27(0.20) 


-0.02 


Jo. 02) 


36.76(0.16) 


1 




1992bo 


0.0172 


15.75(0.13) 


-2.68(0.18) 


0.03 


0.02) 


33.93(0.20) 


1 




1992bl 


0.0422 


17.30(0.08) 


-1.95(0.24) 


0.02 


'0.04) 


35.61(0.17) 


1 




1992bh 


0.0453 


17.58(0.05) 


-0.02(0.25) 


n 1 n 




35.91(0.17) 


1 




1992bg 


0.0365 


16.71(0.09) 


-0.66(0.22) 


0.01 


Jo. 03) 


35.18(0.17) 


1 




1992bc 


0.0196 


15.07(0.11) 


0.51(0.10) 


-0.05 


0.01) 


33.86(0.18) 


1 




1992aq 


0.1009 


19.27(0.05) 


-1.36(0.43) 






37.73(0.17) 


1 




1992al 


0.0135 


14.44(0.16) 


-0.41(0.13) 


-0.05 


0.01) 


33.10(0.22) 


1 




1992ag 


0.0273 


16.26(0.09) 


-0.13(0.21) 


0.19 


0.02) 


34.35(0.18) 


1 




1992ae 


0.0746 


18.39(0.04) 


-0.68(0.21) 


0.01 


0.03) 


36.86(0.16) 


1 




1992p 


0.0265 


16.03(0.09) 


0.70(0.58) 


-0.02 


'0.02) 


34.75(0.19) 


1 




1990af 


0.0499 


17.75(0.05) 


-2.65(0.26) 


0.07 


'0.02) 


35.84(0.16) 


1 




1990o 


0.0306 


16.19(0.08) 


0.39(0.29) 


-0.00 


Jo. 03) 


34.82(0.17) 


1 




1992bk 


0.0589 


18.04(0.13) 


-2.11(0.32) 


-0.04 


'0.08) 


36.35(0.17) 


1 


P 


1992au 


0.0603 


18.12(0.11) 


-1.80(0.53) 


0.01 


Jo. 09) 


36.33(0.17) 


1 


P 


1992j 


0.0461 


17.70(0.09) 


-1.86(0.27) 


0.15 


'0.05) 


35.55(0.16) 


1 


P 


1991ag 


0.0139 


14.37(0.17) 


0.77(0.22) 


0.02 


JO.03) 


32.88(0.22) 


1 


p,z 


1991u 


0.0324 


16.40(0.09) 


0.23(0.22) 


0.11 


JO.03) 


34.62(0.17) 


1 


p 


1991s 


0.0561 


17.71(0.07) 


0.24(0.23) 


0.02 


JO.03) 


36.17(0.15) 


1 


p 


1990y 


0.0387 


17.38(0.08) 


-0.33(0.26) 


0.25 


JO.03) 


35.17(0.17) 


1 


p 


1990t 


0.0397 


17.18(0.07) 


-0.59(0.16) 


0.08 


JO.03) 


35.38(0.16) 


1 


p 


2001CZ 


0.0163 


15.02(0.14) 


0.04(0.14) 


0.12 


JO.Ol) 


33.31(0.14) 


2 




2001cn 


0.0154 


15.26(0.14) 


-0.79(0.12) 


0.22 


JO.Ol) 


33.20(0.15) 


2 




2001bt 


0.0144 


15.26(0.15) 


-1.14(0.11) 


0.26 


JO.Ol) 


33.06(0.16) 


2 


z 


2001ba 


0.0305 


16.19(0.07) 


0.00(0.14) 


-0.04 


JO.Ol) 


34.87(0.09) 


2 




2000ca 


0.0245 


15.53(0.09) 


0.35(0.20) 


-0.07 


JO.Ol) 


34.33(0.10) 


2 




2000bh 


0.0240 


15.92(0.10) 


-0.55(0.17) 


0.08 


JO.03) 


34.24(0.11) 


2 




1999gp 


0.0260 


16.00(0.09) 


1.64(0.21) 


0.06 


JO. 02) 


34.63(0.11) 


2 




1999dk 


0.0139 


14.81(0.16) 


0.29(0.27) 


0.11 


JO. 02) 


33.14(0.17) 


2 


z 


1999cp 


0.0104 


13.93(0.21) 


-0.05(0.25) 


-0.01 


JO. 02) 


32.54(0.22) 


2 


z 


1999cl 


0.0087 


14.81(0.25) 


-0.36(0.20) 


1.20 


JO.Ol) 


30.35(0.26) 


2 


z 


2000ce 


0.0165 


17.02(0.14) 


0.65(0.24) 


0.61 


JO.03) 


34.03(0.15) 


2 


p 



Note. — This tabic is published in its entirety in the eleetronic edition of the Astro- 
physical Journal. A portion is shown here for guidance regarding its form and content. 

P first phase greater than six days after B-band maximum 
^ redshift < 0.015 
° outlier 

number of observations < 5 
' bad light curve fit 

1 Hamuy et al. (1996), ^ Krisciunas et al. (2005), ^ Riess ot al. (1999), Jha ct al. (2006), 
^ Kowalski et al. (2008), Hicken et al. (2009a), ^ Holtzman et al. (2008), ® Schmidt 
et al. (1998, + HZT), " Perlmutter et al. (1999), ^° Barris et al. (2004), " AmanuUah 
et al. (2008), Knop et al. (2003), " Astier et al. (2006), " Miknaitis et al. (2007), 
1^ Tonry et al. (2003), Riess et al. (2007), This Paper 

In this way, all of the 9^ do not have to be explicitly included in the x^i ^ long as the weight matrix, U~^, is 
updated appropriately. Note that may depend on 6^ . 

Minimization over xf^'^ and c*'""'' 

When errors in the independent variable are present, the true values must be solved for as part of the fit (see 
discussion in K08). For one supernova, y = {niB, xi,c), while the model is given by 

{Mb + cosmology) - {ax\''''' - ,5c*™'=), a;f c'™'=) . (C7) 

Using 0^ = (x^i^'^ , c^^'^^) {Mb, a and /3 are global parameters and cannot be handled one supernova at a time), we 
have 




H={ 1 . (C8) 



Spectra and Light Curves of Six SNe la and Union2 33 

The new for this specific example (equation (C4)) is given by equations (4) and (5). 

Minimization over Systematic Errors 

Suppose that we have two kinds of measurements: supernova distance measurements {N supernovae), and zeropoint 
measurements (M zcropoints). The triic valiic of each zeropoint is given by ZPtruo = ZPobserved + AZP. Because the 
uncertainties on the supernova distances are unrelated to the uncertainties on the zeropoints, 

y=( 



Vazp 

where is an x iV block and Vazp is an M x M block. Similarly, y can be split into distances and zeropoints 
y — {Vix^Vazv = 0)- The model for each supernova distance is given by 

Mb + li{z, cosmology) - ^ gzP — ^^^Pa . (C9) 



A 

It is derived in the same way as in the previous section, with the substitutions xi — >■ Xi + ^^^AZP^ and 
c — c + g|p^ AZP A- The model for the true offset of each zeropoint is simply AZP. Thus, 



The upper block oiU ^ is given by 

f^u"pU block = KT' - ^m"'^p(^azp + HlV-'H,)-'HlV-' ■ (CIO) 

none of the other blocks enter into the x^. Inverting this block (see Hager (1989)) gives: 

U = V^,+H^V^-z^Hl (Cll) 

Equation (Cll) gives us the terms we must add to the supernova distance modulus errors to correctly take the 
zeropoints into account. 



